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1 Introduction: Background and Motivations 

Quantum field theory (QFT) is widely viewed as one of the most successful theories 
in science — it has predicted phenomena before they were observed in naturej, and 



£D \ its predictions are believed to be confirmed by experiments to within an extraor- 

dinary degree of accurac}Q. Though it has undergone a long and complex devel- 
opment from its origingj in the 1929/30 papers of Heisenberg and Pauli [105(1106] 
and has attained an ever increasing theoretical sophistication, it is still not clear 
in which sense the physically central quantum field theories such as quantum elec- 
trodynamics (QED), quantum chromodynamics (QCD) and the Standard Model 
(SM) are mathematically well defined theories based upon fundamental physical 
principles that go beyond the merely ad hoc. Needless to say, there are many 
physicists working with quantum field theories for whom the question is of little 
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lr This is an expanded version of an article commissioned for UNESCO's Encyclopedia of Life 
Support Systems (EOLSS). 

2 For example, the existence and properties of the W and Z bosons, as well as the top and 
charm quarks, were predicted before they were found experimentally. 

3 For example, the two parts in one billion difference between the theoretical prediction from 
the Standard Model and the experimentally measured value of the anomalous magnetic moment 
of the muon [H] . 

4 A certain amount of arbitrariness and personal taste must go into pointing to a single point 
of origin, since the 1927 discussion of a quantum theory of electromagnetic radiation by Dirac 
as well as the studies of relativistic wave mechanics by Dirac, Schrodinger and even de Broglie 
were influential. In any case, the interested reader should see |169j for a detailed account of the 
birth of QFT. 
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to no interest. But there are also many who are not satisfied with the concep- 
tual/mathematical state of quantum field theory and have dedicated entire careers 
to an attempt to attain some clarity in the matter. 

This is not the place to explain the grounds for this dissatisfaction; instead, 
the goal of this paper is to provide a perspective on "constructive quantum field 
theory" (CQFT), the subfield of mathematical physics concerned with establishing 
the existence of concrete models of relativistic quantum field theory in a very 
precise mathematical sense and then studying their properties from the point of 
view of both mathematics and physics. Although the insights and techniques 
won by the constructive quantum field theorists have proven to be useful also in 
statistical mechanics and many body physics, these further successes of CQFT are 
not discussed here. In addition, we shall restrict our attention solely to relativistic 
QFT on d dimensional Minkowski space, d > 2; to this point in time most work in 
CQFT has been carried out precisely in that context. Throughout, as is customary 
in QFT, we adopt physical units in which c = h/27f = 1. 

In the 1950's and early 1960's various "axiomatizations" of QFT were formu- 
lated. These can be seen to have two primary goals — (1) to abstract from heuristic 
QFT the fundamental principles of QFT and to formulate them in a mathemati- 
cally precise framework; (2) on the basis of this framework, to formulate and solve 
conceptual and mathematical problems of heuristic QFT in a mathematically rig- 
orous manner. As it turned out, the study and further development of these axiom 
systems led to unanticipated conceptual and physical breakthroughs and insights, 
but these are also not our topic here. 

The first and most narrow axiomatization scheme of the two briefly discussed 
here is constituted by the Wightman axioms (see e.g. |179j ). This axiom system 
adheres most closely to heuristic QFT in that the basic objects are local, covari- 
ant fields acting on a fixed Hilbert space. A (scalar Bose) Wightman theory is 
a quadruple (0, H, U, Q) consisting of a Hilbert space "H, a strongly continuous 
unitary representation U of the (covering group of the) identity component V\ of 
the Poincare group acting upon "H, a unit vector Q £ H which spans the subspace 
of all vectors in T-L left invariant by £7('P+)£| and an (unbounded) operator val- 
ued distribution.!! such that for every test function /, the operator <fi(f) has a 
dense invariant domain T> spanned by all products of field operators applied to Q. 
These conditions are a rigorous formulation of tacit assumptions made in nearly 
all heuristic field theories. In addition, a number of fundamental principles were 
identified and formulated in this framework. 

Relativistic Covariance: For every Poincare element (A, a) £ V\ one has 

U(A, a)4>(x)U(A, a) -1 = <p(Ax + a), in the sense of operator valued distributions 
on T>. 



5 This condition, referred to as the "uniqueness of the vacuum," is posited for convenience. 
With known techniques one can decompose a given model into submodels that satisfy this con- 
dition as well as the remaining conditions 7,22,50,51 . 

6 Although it is possible, indeed sometimes necessary, to choose other test function spaces, 
here we shall restrict our attention to the Schwartz tempered test function space S(R rf ). 



Einstein Causality!] For all spacelike separated x,y G M 4 one has <f>(x)<f>(y) = 
(f>(y)(f>(x) in the sense of operator valued distributions on V. 

The Spectrum Condition (stability of the field system): Restricting one's at- 
tention to the translation subgroup IR 4 C V\, the spectrum of the self-adjoint 
generators of the group f/(R 4 ) is contained in the closed forward lightcone V + = 
{p= (po,pi,P2,Pa) e M 4 | Po-Pi ~vl ~V\ > 0}. 

The reader is referred to [120N179] for a discussion of the physical interpretation 
and motivation of these conditions. There is an equivalent formulation of these 
conditions in terms of the Wightman functions [179] 

W n (x 1 ,x 2 , ...,x n ) = (fi, (j){x 1 )(j){x2) ■ ■ ■ <p(x n )tt) , n e N , 

which are distributions on <S(IR dn ). These two sets of conditions are referred to 
collectively as the Wightman axioms. There are closely related sets of conditions 
for Fermi fields and higher spin Bose fields |120[ll79j . 

A more general axiom system which is conceptually closer to the actual oper- 
ational circumstances of a theory tested by laboratory experiments is constituted 
by the Haag-Araki-Kastler axioms (HAK axioms), also referred to as local quan- 
tum physics or algebraic quantum field theory (AQFT). Although more general 
formulations of AQFT are available, for the purposes of this paper it will suffice 
to limit our attention to a quadruplcl ({A(0)}oeiii T~L, U, Vi) with "H, U and fi as 
above and {A(0)}oen a ne t °f von Neumann algebras A{0) acting on "H, where 
O ranges through a suitable set 1Z of nonempty open subsets of Minkowski space. 
The algebra A{0) is interpreted as the algebra generated by all (bounded) observ- 
ables measurable in the spacetime region O, so the net {A(0)}oeii is naturally 
assumed to satisfy isotony: if 0\ C O2, then one must have A(Oi) C A{C>2)- In 
this framework the basic principles are formulated as follows. 

Relativistic Covariance: For every Poincare element (A, a) G V + and spacetime 
region O G 11 one has U(A, a)A(0)U(A, a)~ l = A(AO + a). 

Einstein Causality o For all spacelike separated regions 0\,Oi G 1Z one has 
AB = BA for all A G A(O x ) and all B G A(0 2 ). 

The Spectrum Condition (stability of the field system): Same as above. 

The reader is referred to [91 I104J for a discussion of the physical interpretation and 
motivation of these conditions. The relation between the Wightman axioms and 
AQFT is well understood (see e.g. [23 l lM l l52.69.183j). It is important to note that, 
in general, infinitely many different fields in the sense of the Wightman axioms 
are associated with the same net of observable algebras. Indeed, an analogy has 



7 Also called microscopic causality, local commutativity or, somewhat misleadingly, locality. 

8 In point of fact, these conditions actually describe an algebraic QFT in a (Minkowski space) 
vacuum representation. By no means is AQFT limited to such circumstances; some other repre- 
sentations of physical interest are briefly discussed below. Moreover, the algebraic approach to 
QFT has proven to be particularly fruitful in addressing conceptual and mathematical problems 
concerning quantum fields on curved spacetimes. 

9 Also often referred to as locality. 



often been drawn between the choice of a particular coordinate system, made in 
order to carry out a computation more conveniently, in differential geometry and 
the choice of a particular field out of the many fields associated with a given net. 
For this and other reasons, those who work in mathematical QFT consider nets of 
observable algebras to be more intrinsic than the associated quantum fields, which 
are used primarily for computational convenience. 

Associated to any Wightman system (0, H,U, fl) is a net of *-algebras V(0), 
O C M 4 . Because all field operators have the common, dense domain T> C 1-t, 
arbitrary "polynomials" of field operators can be formed on T>. V(0) denotes the 
algebra formed by all polynomials (in the sense of functions of infinitely many 
variables) in which the supports of all test functions of all field operators entering 
into the polynomial are contained in the spacetime region O. The algebras V(0) 
are not C*-algebras but satisfy all of the other HAK axioms. Despite the non- 
intrinsic nature of such algebras and despite the technical disadvantages of working 
with *-algebras instead of with C*-algebras, mathematical quantum field theorists 
find it convenient for various purposes to work with such nets or with similar nets 
of non-C*-algebras. 

The goal of constructive QFT, as is it usually understood, is to construct in a 
mathematically rigorous manner physically relevant quantum field models which 
satisfy one of these systems of axioms and then to study their mathematical prop- 
erties with an emphasis on those properties which can be shown to have physical 
relevance. This article briefly describes such models and the means by which 
they were constructed and is organized both historically and by the construction 
techniques employed. 

As pointed out independently by Borchers and Uhlmann, the Wightman axioms 
can be understood in a representation independent manner in terms of what is now 
called the Borchers (or Borchers-Uhlmann) algebra — a tensor algebra constructed 
out of the test function space 5(IR 4 ) with operations directly motivated by the 
Wightman axioms. Borchers algebras have been extensively studied from the 
point of view of QFT, especially by Borchers, Uhlmann, Yngvason and Lassner 
(see e.g. [109] for definitions and references). A Wightman system can be thought 
of as a concrete representation of the Borchers algebra, and for a time there was 
hope one could arrive at quantum field models by defining suitable states on the 
Borchers algebra and employing the standard GNS construction to obtain the 
corresponding representation. However, it proved to be too difficult to conjure 
such states. 

The first quantum field models constructed were the free quantum fields, the 
Wick powers of such free fields and the so-called generalized free fields. These 
models have been constructed using a variety of techniques (cf. e.g. [E1II21EH1ISI2 
I120yi45pl74yi96j ) and have been shown to satisfy the two axiom systems discussed 
above; a recent construction of free fields which is of particular conceptual interest 
is briefly described in Section |6j The Hilbert space upon which such fields act is 
called the Fock space. Common to these models is the fact that their S-matrix, 



the object which describes the scattering behavior of the "particles" described by 
such fields (cf. [9"| I12U] ). is just the identity map. 

We turn now to models with nontrivial S-matrices, i.e. interacting quantum 
field models. When referring to the models, we employ the standard notation 
Md, which means quantum model M in d spacetime dimensions. Because the 
mathematical and conceptual difficulties inherent in the construction of quantum 
field models are quite daunting, constructive quantum field theorists proceeded by 
considering increasingly challenging models; this often entailed starting the study 
of the model M with d = 2, then d = 3, and finally d = 4. At this point in time 
only a few models have been constructed in four spacetime dimensions. In this 
respect, the reader is referred to Section [S] for a few words about the outlook for 
CQFT after nearly fifty years of strenuous effort. The reader should note that 
all results discussed in this article, unless explicitly stated otherwise, are proven 
according to the criteria accepted by mathematicians and not merely on the basis 
of the plausibility arguments accepted by most physicists as "proof" . 



2 Algebraic Constructions I 

Preceded by the 1965 dissertations of Jaffe |116j and Lanford [132] . the first con- 
structions of interacting quantum fields were carried out in the late 1960's and 
early 1970's. In this early work the real time models were constructed directly us- 
ing operator algebras and functional analysis as the primary tools. Due to Haag's 
Theorem, it was known that the Hilbert space in which these interacting quantum 
fields would be defined could not be Fock space (see e.g. P-84J). However, because 
Fock space was the sole available starting point at that time, "cutoff's" were placed 
on the interacting theories so that they could be realized on Fock space in a math- 
ematically meaningful manner. These cutoffs were of two general kinds -- finite 
volume cutoffs and ultraviolet cutoffs -- each addressing independent sources of 
the divergences known in QFT since early in its development. Guided by heuris- 
tic QFT's division of Lagrangian quantum field models into superrenormalizable, 
renormalizable and nonrenormalizable modelsjlj the constructive quantum field 
theorists began with the simplest category, the superrenormalizable models. To 
be able to address the infinite volume divergence without wrestling simultaneously 
with the ultraviolet divergence, constructive quantum field theorists first consid- 
ered self-interacting bosonic quantum field models in two spacetime dimensions. 

We begin with Glimm and Jaffe's construction of the (</> 4 )2 model [STH54] . the 
self-interacting scalar Bose field on two dimensional Minkowski space with La- 
grangian self- interaction A0 4 , where A is the coupling constant. Let "Ho be the Fock 
space for a (free) scalar hermitian Bose field (f)(t, x) of mass m > ((£, x) G 1R 2 ). Let 
7r(£, x) = d(p(t, x)/dt be the canonically conjugate momentum field and V C "Ho be 
the dense set of finite-particle vectors in "H - Then, for every / in a dense subspace 



ln This classification is based upon the perturbation theory associated by Feynman and others 
with interacting fields, viewed as perturbations of free fields. 
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of I/ 2 (M), the operator (f>o(f) = f (p(0,x)f(x) dx is essentially self-adjoint on 
T> and (f) (f)V C T> (similarly for 7r (/)). These operators satisfy the canonical 
commutation relations (CCR) on T>: 

Mf)Mg) - Mg)Mf) = i<f,g>*, 

Mf)M9) - M9)Mf) = = n (f)n (g) - 7r (g)7r Q (f) , 

for all f,g G «S(M), where < •, • > is the inner product on L 2 (M) and 1 is the 
identity operator on %. When exponentiated using the spectral calculus, (the 
closures of) these operators provide a Weyl representation of the CCR. For each 
bounded open subset Oct, denote by .4(0) the von Neumann algebra generated 
by the Weyl unitaries 

{e^ o(/) ,e^ o{/) | / G <S(R) , supp(/) C O} . 

(supp(/) denotes the support of the function /.) 
The total energy 

H = - : (tt(0, x) 2 + V0(O, x) 2 + m 2 0(O, x) 2 ) : dx 

of the free field is a positive quadratic form onDxD and therefore determines 
uniquely a positive self-adjoint operator, which we also denote by H . The double 
colons indicate that the expression between them is Wick ordered, which is a 
physically motivated way to define in a rigorous manner a product of operator 
valued distributions. In this case, the Wick ordering is performed with respect 
to the Fock vacuum (c/. [SD])- With g G L 2 (M) nonnegative of compact support, 
Glimm and Jaffe showed that, for each A > 0, the cut-off interacting Hamilton 
operator 

H(g) = H + A J : 0(0, x) 4 : g(x) dx 

is essentially self-adjoint on T>\ n \ and its self-adjoint closure, also denoted by H(g), 
is bounded from below. By adding a suitable multiple of the identity we may take 
to be the minimum of its spectrum. Then, they proved that is a simple eigenvalue 
of H(g) with normalized eigenvector fi(g) E "Ho- 

For any { 6 R, let O t denote the subset of M consisting of all points with 
distance less than \t\ to O. By choosing the cutoff function g to be equal to 1 on 
Of, then for any A G 4.(0) the operator 

a t (A) = e ltH ^Ae- itH{9) 

is independent of g and is contained in 4(O t ). For any bounded open 0cl 2 and 
t G R, let Q(t) = {x G R | (t,x) G O} be the time t slice of O. We define A{0) 



11 Without the cutoff g, the interacting Hamilton operator is not densely defined in Fock space. 



to be the von Neumann algebra generated by [j s <7 s (*4(0(s)))o Finally, we let A 
denote the closure in the operator norm of the union [J A(O) over all open bounded 
(9cR 2 . Hence, at is an automorphism on A and implements the time evolution 
associated with the interacting field. Similarly, "locally correct" generators for 
the Lorentz boosts and the spatial translations can be defined, resulting in an 
automorphic action a on A of the entire (identity component of the) Poincare 
group V\ in two spacetime dimensions. 

For each A G A, we set uj 9 (A) = < Q(g), AQ(g) > ((•, •) denotes here the inner 
product on %) to define the locally correct vacuum state u g of the interacting field. 
Taking a limit as the cutoff function g approaches the constant function 1, Glimm 
and Jaffe showed that uj g (A) — > u(A), for each A G A, defines a new (locally 
normal) state u on A which is Poincare invariant, i.e. u(a^ x ^(A)) = u>(A) for all 
(A, x) G V\ and all A G A. Employing the GNS construction, one then obtains a 
new Hilbert space 7i, a representation p of A as a C*-algebra acting on H, and a 
vector Q G 7i such that p(A)Q is dense in 7i and 

uj(A) =< Q, p(A)Q > , for all A G A . 

In addition, one obtains a strongly continuous unitary representation U of the 
Poincare group in two spacetime dimensions under which the algebras p(A{0)) 
transform covariantly, i.e. 

[/((A, x)) p(A(0)) C/((A, x))- 1 = p(A(AO + x)) . 

Both the HAK and Wightman axioms have been verified for this model. 

The generators of the strongly continuous Abelian unitary groups 
{p( e ^(/)) | t G R} and {p(e it7r(/) ) | t G R} satisfy the CCR. However, this repre- 
sentation of the CCR in % is not unitarily equivalent to the initial representation in 
Fock space, in accordance with Haag's Theorem. Indeed, by taking different values 
of the coupling constant A in the above construction, one obtains an uncountably 
infinite family of mutually inequivalent representations of the CCR (see [70J). 

It is in this representation (p, %) that the field equations for this model find a 
mathematically satisfactory interpretation, as shown by Schrader |166j . And it is 
to the physically significant quantities in this representation that the corresponding 
perturbation series in A is asymptotic — see below for further discussion. For this 
and other reasons, u is interpreted as the exact vacuum state in the interacting 
theory corresponding to the Lagrangian interaction A0 4 , and the folium of states 
associated with this representation contains the physically admissible states of the 
interacting theory. Many further properties of physical relevance have been proven 
for this model more recently — see the discussion below in Section 13.11 

The results attained for the 0^ model were subsequently extended to -P(0)2 
models (using a periodic box cutoff), where P{4>) is any polynomial bounded from 



12 One can then show that the algebra A(0) coincides with the von Neumann algebra generated 
by bounded functions of the self-adjoint field operators J <p(t,x)f(t,x)dxdt, with test functions 
f(t,x) having support in O. 



below [8"oT - [8"8"] ] 13 l (See [8"M?U] for more complete references and history of this devel- 
opment.) Hoegh-Krohn employed the techniques of Glimm and Jaffe to construct 
models in two spacetime dimensions (with similar results) in which the polyno- 
mial interaction P(4>) is replaced by a function of exponential type, the simplest 
example being e Q< ^ [107j . 

Analogous results were proven for Y 2 , the Yukawa model in two spacetime 
dimensions by Glimm and Jaffe and Schrader [85|l89l[T65] . In this model one 
commences with the direct product T-Lq = Tib <H> %f of the Fock space Tib for a 
scalar hermitian Bose field <fi(t, x) of mass nib > and the Fock space Tif for a 
Fermi field ip(t, x) of mass m/ > 0. In this model the free Hamiltonian H is the 
total energy operator of the free fields (f) and ip. Because there is still an ultraviolet 
divergence remaining after Wick ordering, the cutoff interacting Hamiltonian is 
H (g, k) = H + Hj(g, K) + c(g, k), where Hj(g, k) is the result of applying a certain 
multiplicative ultraviolet cutoff (which is removed in the limit k — > oo) to the 
formal expression 

Hi{g) = A / g(x) 0(0, x) : iptp : (0, x) dx , 

and c(g, k) is a (finite) renormalization counterterm determined by second-order 
perturbation theory which diverges as k — > oo [85] and includes both a mass and 
vacuum energy renormalization. With both volume and ultraviolet cutoffs in place, 
H(g,n) is a well defined operator on 1-Lq. Glimm and Jaffe show that as k — > oo the 
operator H(g, k) converges in the sense of graphs to a positive self-adjoint operator 
H(g) with an eigenvector Q(g) G T-Lq of lowest energy 0. Once again, they control 
the limit as g — > 1 of the expectations o-> g {A) for all A G Ab® Af and obtain a 
state u on Ab <8> Af that provides a corresponding (GNS) representation of the 
fully interacting theory. Glimm and Jaffe [89J also proved that the Yukawa field 
equations are satisfied by the fields in that representation. A similar argument 
was applied to the Y 2 + -P(0)2 model by Schrader |164j ). where 

Hj{g) = A [ g{x) (0(0, x) : ^ : (0, x) + : P(</>) : (0, x)) dx 

and P{4>) is any polynomial bounded from below. The axioms of HAK and Wight- 
man were shown to hold in these models (see p22]), at least for all sufficiently small 
values of the coupling constant A. In addition, by using a mixture of algebraic and 
Euclidean methods Summers |182] showed that the model manifests further prop- 
erties of physical relevance, such as the funnel property (also known as the split 
property) and all assumptions of the Doplicher-Haag-Roberts superselection the- 
ory (c/. [9 tll04j ). Therefore the model also admits the physically expected Poincare 
covariant, positive energy, charged representations associated with the generator 
of the global U(l) gauge group of the model, which are mutually unitarily inequiv- 
alent. 



13 If P((f) is not bounded from below, then the corresponding cutoff Hamiltonian H (g) is not 
bounded from below and the resulting model is not stable. 
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Along the lines employed in the construction of the Yukawa model in two space- 
time dimensions, Glimm and Jaffe [91] also showed for the 0| model that the spa- 
tially cutoff Hamiltonian H(g) is densely defined, symmetric and bounded below 
by a constant E(g) proportional to the volume of the support of g. The renormal- 
ization constants in the Hamiltonian H (g) are again given by perturbation theory 
and involve counterterms to the vacuum energy and the rest mass of a single par- 
ticle. The proof was technically more challenging than that for Y2, even though 
the results were more limited to a significant extent. There was real motivation to 
find an alternative approach, as described in the next section. 

However, before proceeding to the next section we mention the Federbush 
model, a model of self-interacting fermions in two spacetime dimensions. First 
proposed by Federbush [61], the Lagrangian of the model is 

^2 4> s ($ ~ m(s))4>s - 2'k\c i1u J^J v _ 1 , 

s=±l 

where e w = — e i = 1, eoo = e n = 0, J^ = ip^ips and m(s) > 0, s = ±1. 
Without cutoffs of any kind, a concrete realization of the Federbush model can 
be given in terms of certain exponential expressions on a suitable Fock space, and 
Ruijsenaars |159] proved that this realization satisfies the Wightman axioms when 
A G (— |, |) (Einstein causality is actually verified only for sufficiently small A). Of 
particular interest, he proved that the associated Haag-Ruelle scattering theory 
is asymptotically complete. The S-matrix is nontrivial, but there is no particle 
production [160] , The Federbush model was the first non-superrenormalizable 
model for which any of these properties have been proven. 



3 Functional Integral Constructions — Euclidean 

The technical difficulties of the approach described in the preceding section were 
formidable, and in the early 1970's a technically more manageable program re- 
placed (at least temporarily) the earlier approach] 14 ! In this newer program a 
(scalar Bose) Euclidean quantum field is a random variable valued distribution 
<p(x) with probability measure fi on the corresponding distribution space. The 
moments of this measure 

S n (xx,x 2 ,...,x n )= / (f)(xi)4>(x 2 ) ■ ■ ■ </)(x n ) dfi((p) , n G N, 



are called the Schwinger functions. Formally, the connection between the Wight- 
man functions and the Schwinger functions is described as follows: Let x = 
(x , xi, . . . , Xd-i) be the coordinates of a point in d dimensional Minkowski space 
and consider x E = (—ix , Xi, . . . , x^-i)- Then 

S n (x!, x 2 , ■ ■ ■ , x n ) = W n (xf, xf , . . . , x® ) , n G N . 



14 See [5511 1741H571I 188] for some history and references. 
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In point of fact, it is known that the Wightman functions of a Wightman theory can 
be analytically continued to such Euclidean points (xf, xf, . . . , x„) |12U|I187| so the 
above relation is meaningful in that sense. Under certain additional conditions the 
resulting Schwinger functions are the moments of a Euclidean invariant probability 
measure /i. For example, the Schwinger functions of the free scalar Bose field of 
mass mo > on two dimensional Minkowski space are the moments of the Gaussian 
probability measure \xc on 5'(1R 2 ) with mean and covariance C = (—A + m 2 ] )~ l . 
The corresponding two-point Schwinger function is 

r „ik-(x-y) 

S 2 (x-y) = (2n)- 2 / ———dk 2 , 



k 2 + wig 

where all inner products are Euclidean. 

Hence, from a Wightman field theory one can obtain in this manner a Euclidean 
field theory. Under certain circumstances one can also obtain a Wightman field the- 
ory from a Euclidean field theory by appropriate analytic continuation of Schwinger 
functions to obtain Wightman functions satisfying the Wightman axioms. Of the 
various sets of conditions known to permit this reconstruction of the real time the- 
ory from the imaginary time theory, we mention only the Osterw •alder- Schrader 
axioms p.52]. Here one commences with a sequence {S n (x\, x 2 , ■ ■ ■ ,x n )} ne ^ of 
tempered distributions (again called Schwinger functions) completely symmetric 
in their arguments which are invariant under the natural action of the Euclidean 
group. One requires additionally of these distributions that they satisfy a cer- 
tain cluster property (ergodicity under the "time" translations), the Osterwalder- 
Schrader positivity condition (also called reflection positivity) and a technical lin- 
ear growth condition. Osterwalder and Schrader showed that an inductive con- 
struction of analytic continuations applied to these Schwinger functions yields a 
set of Wightman functions satisfying the Wightman axioms. They also proved 
that analytically continuing the Wightman functions of a Wightman theory in 
the manner indicated above results in Schwinger functions satisfying all of the 
Osterwalder-Schrader axioms except possibly the linear growth condition. Since 
the linear growth condition may fail, the Osterwalder-Schrader axioms may not 
be equivalent to the Wightman axioms] 15 ! 

In general, a given sequence of Schwinger functions satisfying the Osterwalder- 
Schrader axioms need not be the moments of a measure /i. However, the gen- 
eral strategy of the Euclidean construction program was to construct a suitable 
probability measure /i on S'(R d ) such that its moments satisfy the Osterwalder- 
Schrader axioms. Application of the Osterwalder-Schrader reconstruction theorem 
would then result in a corresponding relativistic quantum field model satisfying 
the Wightman axioms. One would then further study these models for properties 
of physical relevance that go beyond the basic axioms. This approach also stimu- 
lated a fruitful exchange of techniques and tools between quantum field theory on 



15 Finding a physically meaningful set of conditions on Schwinger functions which is equivalent 
to the Wightman axioms is an open problem — see [163) for a recent overview of the topic 
"equivalence of Wightman functions and Schwinger functions" with references. 
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one side and probability theory and classical statistical mechanics on the other. 

The connections between Euclidean QFT and probability theory are even richer 
than what has been suggested above. Nelson |146j pointed out that the free scalar 
massive Bose Euclidean field is a Markov process, as is the finite volume P(4>)2 field. 
This line of thought led to the use of the stochastic processes which are solutions 
of suitable stochastic differential equations to construct models of quantum fields 
— see Section [5j 

3.1 P((f>) 2 models 

The P(<f))2 models were reexamined from the point of view of the Euclidean con- 
struction program, and new results of physical relevance were obtained. In analogy 
to the real time construction sketched in Section [21 one commences with a local 
perturbation of the free field, here represented by the Gaussian measure He- Once 
again, let P(<f>) be a polynomial bounded from below and A be a rectangular re- 
gion in Euclidean M 2 . The Wick ordered powers of the Euclidean free field may 
be defined recursively by the conditions : 0° := 1, d : <p n : jd<\> = n : n_1 : and 
/ : n : dfi C = 0, n G N. Let U A = J A : P(<f>) : (x) d 2 x. Then e~ UA G LP(d^ c ) 
for all p < oo |174j , so 

^ = fe-^dfic ( } 

is a well defined probability measure on <S'(M 2 ), where the coupling constant A > 0. 
One may then consider the corresponding Schwinger functions 



S n (xi, x 2 ,...,x n )= / 0(a?i)0(ar 2 ) • ■ ■ <p(x n ) dfi A , n e N . 

By employing cluster expansions Glimm, Jaffe and Spencer proved that the limits 
S n (x!, x 2 , ■ ■ ■ , x n ) = lim S^(xi, x 2 , . . . , x n ) , n e N 



exist] 16 ! and satisfy the Osterwalder-Schrader axioms for all sufficiently small A/m 2 , 
[92] . Hence, one obtains a corresponding relativistic quantum field model. 

Similar results can also be obtained, at least for a large class of polynomials, by 
commencing with a lattice approximation and controlling the limit as the lattice 
spacing goes to [9511174] and/or by using correlation inequalities instead of cluster 
expansions [95 11174] . In the former, placing the quantum field model on a lattice 
results in a classical Ising ferromagnetic system with unbounded spins and nearest 
neighbor coupling arising from the finite difference approximation to — A + m 2 ,. 
The lattice spacing serves as an ultraviolet cutoff. An example of the latter are 
the GKS-inequalities, which were extended from lattice ferromagnetic systems to 

16 In point of fact, the measures /^a converge weakly to a measure /i m whose moments are 
S n (xi,x 2 ,-.-,x n ), n e N. 
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P(0) 2 models by Guerra, Rosen and Simon |1U2] : For P{4>) = Q(4>) ~ h(f>, where 
Q is an even polynomial and h > 0, one has 

S%{x 1 ,x 2 ,...,x n ) >0,nGN, 

and 

'-'n+ml 3 '!' ^2; • • • > %>n+m) ^ >->« (^1) ^2? • • • i %n)<5 m \%n+li %n+2i ■ ■ ■ i %n+m) > 77, 777 G .N , 

in the sense of distributions. Nelson proved that the Schwinger functions are 
monotone in the volume cutoff Aj^| namely for suitably regular regions A'c A 
one has 

S n (/l, /2, ■ ■ ■ , fn) < S n (/i, /2, • • • , /n) 

for all positive test functions /i,...,/ n and all n G N |146j . Bounds on the 
Schwinger functions which are uniform in A yield the existence of the infinite 
volume limit. The GKS inequalities are used in the proof of both the monotonicity 
and the uniform bounds. Further arguments yield the validity of the Osterwalder- 
Schrader axioms for the limit Schwinger functions. Though use of correlation 
inequalities restricts the interaction polynomial, wherever correlation inequalities 
can be used, there is no restriction on the size of A > 0. This is an advantage with 
respect to the results attained using cluster expansions. 

In the weak coupling limit (sufficiently small A/m 2 ,) the "uniqueness of the 
vacuum" follows from exponential clustering bounds on the Schwinger functions 
(or the Wightman functions): there exists an m > such that for all < k < n G N 

and a G M 2 

S n {Xi, X2, • • • , Xk, Xk+i — a, Xk+2 — a, ■ ■ ■ , x n — a) — (3-2) 

Sk{Xi,X2, ■ ■ ■ ,Xk) S n -.k{Xk+i, x k+2, ■ ■ ■ ,x n ) < Ce m " a " . (3.3) 

In fact, the bound (j3.2p establishes also the existence of a mass gap at least as 
large as m, i.e. the spectrum of the mass operator is contained in {0} U [m, oo). 
Indeed, it was shown that there exists an 777 > so that the spectrum of the 
mass operator is contained in {0, 777} U [mi, 00) and that as A — > with 7% fixed, 
one has m —¥ m and m\ — > 2mo [S2]- Hence, weakly coupled P[^)i models 
have an isolated one-particle hyperboloid, and the Haag-Ruelle scattering theory 
(see [9| ll20j ) may be applied [92], resulting in a well defined scattering theory 
for these models. It has been shown by Osterwalder and Seneor [153j that the 
resultant S-matrix is nontrivial and that the usual perturbation expansion in the 
coupling constant for this S-matrix is asymptotic to the exact S-matrix. In these 
models there is particle production. 

Much work has been expended in the study of the particle structure of P((f>)2 
models which goes beyond the existence of an isolated mass hyperboloid, since this 
is of direct physical interest and also because it appears likely that such knowledge 
would be a necessary prerequisite for any proof of the asymptotic completeness 



17 Strictly speaking, "half-Dirichlet boundary conditions" on A are introduced to prove this. 
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of the models. An early example is the proof by Spencer and Zirilli |177] that 
in the A0 2 model there are no even bound states of energy less than 4(m — e), 
where m = m(A) is the physical mass of the model and e — > as A — > 0. On 
the other hand, Dimock and Eckmann [45] showed that in the A(0 6 — 4 ) 2 model 
there exists a unique two particle bound state for sufficiently small A, and it has 
mass 2m(l — g^jA 2 + 0(A 3 )). They subsequently generalized these results [46] 
to any polynomial interaction having no positive 4 term. They further showed 
that if there is such a 4 term, then there is no two particle bound state. In both 
cases, there are no bound states masses embedded in the continuum below 3m — e, 
i.e. the model is two-particle asymptotically complete. Neves da Silva [147] proved 
that when the interaction polynomial is even and there exists a two particle bound 
state, then there also exists a three particle bound state of energy less than 3m. 
The uniqueness of the former entails the uniqueness of the latter. Further details 
of particle structure are strongly model dependent, but general results attained by 
Bros and Iagolnitzer [27] are a useful starting point. 



Using Euclidean Markov fields and correlation inequalities, Albeverio and Hoegh- 
Krohn [2] proved that the infinite volume limit of scalar bosons with even expo- 
nential self-interaction in two space-time dimensions satisfies the Osterwalder- 
Schrader axioms, including clustering, so that the corresponding Minkowski space 
theory satisfies the Wightman axioms, including uniqueness of the vacuum. To- 
gether with Gallavotti they showed [3] that when d > 3 the model can also be 
constructed, but its Schwinger functions coincide with those of the corresponding 
free field, i.e. it is trivial. 

In general QFT the "uniqueness of the vacuum" can fail in a physically sig- 
nificant manner — the quantum field model can manifest "phase transitions" en- 
tirely analogous to those found in statistical mechanics. In this circumstance the 
subspace Ho of vectors in H invariant under the action of U{V\) is not one- 
dimensional, and one can find two unit vectors fii,fi 2 G "Ho an d an observable 
A (often called the order parameter) such that (Qx,AQi) ^ (fi 2 ,Afi 2 ). As in 
statistical mechanics, one can determine "phase diagrams" plotted against the 
various parameters of interaction in the model with "critical points" where phase 
transition curves end. 

The first such results were proven for P(4>) = (<f) 2 — a 2 ) 2 /a 2 for a a sufficiently 
large constant. By a suitable change of variables and Wick ordering, this is phys- 
ically equivalent to the choice P((f>) = A0 4 + \m 2 ] (f> 2 for sufficiently small mp/Acj 
As in statistical mechanics, one displays the phase transition by considering the 
interaction P(4>) = {(f) 2 — o 2 ) 2 ja 2 — h(j), with h a constant representing an "external 
field" . For h ^ the corresponding model has a unique vacuum, by the Lee- Yang 
theorem generalized to the 2 m odel by Simon and Griffiths [173] , i.e. the corre- 
sponding measure \x a ^ is ergodic. For all sufficiently large a, Glimm, Jaffe and 



18 In direct analogy to statistical mechanics, the limit rn^/X — > is referred to as the low 
temperature limit and that of A/toq — > is called the high temperature limit. 
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Spencer show that 



lim (<j>(f)) a , h = - lim (<j>(f)) a , h > 

h\0 h/-0 



for positive test function /, where (-)a,h is the expectation with respect to n a> h- 
Hence, there is a phase transition at h = for sufficiently large a and the field 
itself is an order parameter. Indeed, Glimm, Jaffe and Spencer [91] subsequently 
modified the cluster expansion to show that with P(0) = A0 4 — ^<fi 2 — hep — E c (E c 
chosen so that the infimum of the polynomial is 0) and \x\, the Gaussian measure 
on <S'(IR 2 ) with mean b = ±(8A)~ 1//2 (note these are the minima of P(<j))) and 
covariance C = (—A + 1) _1 , then if h and b have the same sign the Schwinger 
functions associated with the measure 

e - f A (:P(0:(*)+£:(*-6) a :(*))rf»a! ^ 

^ A = J e ~ / A (■■P(<t>y(x) + \:(<t>-b?:(x))d^x dflb { ■ > 

converge as A /* IR 2 to Schwinger functions satisfying the Osterwalder-Schrader 
axioms for all sufficiently small A, in particular, the corresponding vacuum is 
unique and the limit measure y,± is ergodic. Note that in (j3.4p the mean and mass 
in the Gaussian measure exactly cancel the term added to P(4>) in the region A, 
leaving a ± boundary condition outside of A. What is more, (<f>(x))± = &+0(A 3 / 2 ), 
where (-)± is the expectation in fi±. Hence, for h = one obtains two distinct 
Wightman theories in which the corresponding vacuum expectations of the field 
differ by a sign and in which the <fi — > —<p symmetry of the action is spontaneously 
broken. In analogy to situations arising in statistical mechanics, one says that at 
h = there are two distinct "phases" for the model. Koch |129] showed that in 
each of these pure phases the mass hyperboloid is isolated, so that Haag-Ruelle 
scattering theory may be applied, and that there exists a two particle bound state. 

The mean field cluster expansion of Glimm, Jaffe and Spencer was subsequently 
extended by Summers [181] to an interaction involving a sixth degree polynomial, 
which resulted in a value of the interaction coefficients at which three distinct 
phases coexist and in the corresponding phase diagram there are phase lines which 
are not associated with spontaneous symmetry breaking. There are also values 
of the coefficients resulting in two "critical points" where phase lines simply end, 
i.e. where the distinction between the phases vanishes. Each of the pure phases 
satisfies all of the Osterwalder-Schrader axioms, and the perturbation series for 
the Schwinger functions (suitably adjusted about the mean value of the field in 
the phase) is asymptotic to the exact (similarly adjusted) Schwinger functions. 
Imbrie [113] generalized this study of phase transitions in P(<p)2 models to very 
general polynomials yielding quite complex phase diagrams by suitably adapting 
techniques developed in statistical mechanics to this situation. 

Although most effort has been expended on the construction of quantum field 
models in a vacuum representation, other representations are of physical interest, 
as well. One such class of representations studied by constructive and axiomatic 
quantum field theorists is the class of equilibrium thermal representations, i.e. 
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representations associated with states which satisfy the KMS condition for some 
fixed inverse temperature (cf. [26J for further background and references). Here 
we only mention results proven in concrete interacting models. Hoegh-Krohn [11)8] 
considered scalar Bose fields in two spacetime dimensions with either polynomial or 
exponential self-interaction and showed using a hybrid of algebraic and Euclidean 
techniques that for any T > the infinite volume limit of the spacetime cutoff 
Gibbs state at temperature T exists on a suitable global algebra of observables 
and satisfies both the KMS condition for (3 = T _1 and the cluster property and is 
translation invariant. More recently, Gerard and Jakel [77] revisited the matter and 
provided a different proof of these results, as well as some further, more technical 
observations. 

Thermal equilibrium states cannot be Lorentz invariant [150] , even if the equa- 
tions of motion are invariant under Poincare transformations and the signal prop- 
agation speed is finite. Nonetheless, as recognized by Bros and Buchholz [28] . 
the passivity (cf. [26]) of an equilibrium state should still be visible to an ob- 
server in motion with respect to the rest frame distinguished by the KMS state. 
They showed that this entails that the thermal equilibrium states of a relativis- 
tic QFT should have stronger analyticity properties in configuration space than 
those already implicit in the KMS condition. They formulated these properties 
as a relativistic KMS condition. Using the Euclidean approach to thermal fields 
and technical advances in treating spatially cutoff models due to Klein and Lan- 
dau [127] . Gerard and Jakel [78] proved that the two-point function in the P{4>)2 
model satisfies this relativistic KMS condition, and Jakel and Robl [115] extended 
the result to general n-point functions in this model. 

Other non- vacuum representations of physical interest which have been studied 
rigorously in concrete models are charged representations, an example of which 
is provided by the Yukawa2 model [182] . where the countably infinite eigenspaces 
of a global "charge" operator Q commuting with all observables and the repre- 
sentation of the Poincare group yield mutually inequivalent representations of the 
observables ( "superselection sectors" ) that satisfy the HAK axioms, excepting the 
existence of the vacuum (the vacuum sector is the charge eigenspace of Q). 

Soliton representations can arise with the appearance of phase transitions and 
associated spontaneous breaking of symmetries of the model, as in the $\ models 
discussed above. They are believed to occur primarily in models in two spacetime 
dimensions, though there are heuristic indications that they might exist in non- 
abelian Yang-Mills models in four spacetime dimensions. They are Poincare co- 
variant, positive energy representations of the fields (observables) that, in a certain 
sense, interpolate between two distinct vacuum representations. As an illustration, 
consider the <p\ model in the parameter range where the phase transition discussed 
above occurs. At h = there exist two vacuum vectors Q± inducing two vacuum 
states u>±(-) = (fi±, -fl±) on the fields, equivalently on the algebra of observables 
A. As shown by Frohlich [71] . there exists a localized automorphism a on A such 
that u)± o a are states on A and the GNS representations of A corresponding to 
w± o a are Poincare covariant and satisfy the spectrum condition. If r H.± are the 
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representation spaces of A corresponding to u>± and % SjS / the representation spaces 
corresponding to u± o a, then Ti+ © %- © % s © "H s ' is the full Hilbert space of 
the 02 model (in the stated parameter range). On this space a time- independent 
"topological" charge operator is defined formally by 



f d 

Q = —4>(t,x) dx. 



H + ®T-i- is the eigenspace of Q corresponding to the eigenvalue 0, and H s , resp. 
7i s i, is the eigenspace corresponding to the eigenvalue 2(f2 + , 0(x)f2 + ) > 0, resp. 
2($1_, 0(x)fi_) = —2(Q + ,(j)(x)Q + ). The sense in which the soliton, resp. anti- 
soliton, state u s , resp. uv, interpolates between the two vacua is suggested by 



and 



lim ou s ((f)(t,x)) = u>±((j)(t,x)) 

:— >±oo 



lim u s '((j)(t,x)) = u) T ((j)(t, x)) . 

x— >±oo 



Although not all details have been published, there are strong indications in the 
work of Frohlich and Marchetti [7HE1] that there are single soliton (antisoliton) 
states in 7-L s (*H s r) which are created out of the vacuum vectors Q + (fi_) by a local 
field s(x), called the soliton field. Bellisard, Frohlich and Gidas [18] have proven 
that for small enough coupling A the mass of the soliton grows like A -1 . 

These ideas have been explored also in P(<f>)2 models, the Yukawa model and 
others. Of particular interest is the sine-Gordon model (see Section [376]) . which, 
because of the periodicity of the interaction, admits countably infinitely many 
soliton sectors [7TJ. These and other models have also been studied from the 
Euclidean point of view by Frohlich and Marchetti [73], casting a new perspective 
on these results. See also Schlingemann |162j for more recent developments. 



3.2 The <f>\ model 

Also the 0| model was revisited using Euclidean techniques, yielding stronger 
results. As seen in Section [21 the 0g model is superrenormalizable and requires an 
infinite mass renormalization. Hence, both a volume and an ultraviolet cutoff are 
necessary to begin with well defined quantities. Let nc K be the Gaussian measure 
on the Schwartz distribution space iS'(R 3 ) with mean and covariance C K , with 
kernel 

C K (x -y) = (2tt)- 3 / e ' [ \ V {k/K) dk 3 , (3.5) 

mo > 0, where r\ is a smooth function of compact support taking the value 1 in 
a neighborhood of the origin in M 3 . The cutoff free field K is the corresponding 
random variable valued distribution. The cutoff action is 

V{\, A, k) = / A : 0^ : (x) + c{k, A, K (x)) d 3 x , 
Ja 
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where c(k, A, (/> K (x)) is again a renormalization counterterm suggested by pertur- 
bation theory (see e.g. [138] ). and A > is the coupling constant. Define the cutoff 
interacting measure to be 



fi\,A,K = 



e -V(X,A, K ) ^ 

J e-v(AA.) ^ 



Using phase cell and cluster expansions, Feldman and Osterwalder [63J, on one 
hand, and Magnen and Seneor |138j . on the other, have independently shown 
that for sufficiently small \/rn^ the corresponding Schwinger functions converge 
as k — > oo, A /*■ IR 3 , to Schwinger functions satisfying the Osterwalder-Schrader 
axioms, so that there is a corresponding Wightman theory. By using different 
techniques going back to Nelson and Guerra, Seiler and Simon |171] were able to 
remove the restriction on X/ttIq. Using a lattice approximation and correlation 
inequalities, Park |155j attains the same results. The resulting quantum fields 
satisfy the corresponding 4 field equations [64J , are locally associated with a net of 
algebras satisfying the HAK axioms [171] , and the perturbation expansion for the 
Schwinger functions is Borel summable [140] . Burnap [38] showed the existence of 
one particle states in this model, so the Haag-Ruelle scattering theory is applicable 
to this model. 

In the intervening thirty years the model has been revisited from many points of 
view; much effort has been exerted to simplify the techniques used in the original 
proofs (cf. the recent [143] ). and we mention, in particular, that various rigorous 
versions of the renormalization group have been developed and applied to the (f)\ 
model (cf. [32] for an overview). This approach and a second, using self-avoiding 
random walks to construct 0^ models for d > 2, are briefly discussed in Section 



3.3 (\>\ models and their dependence on d 

According to the classification arising from the application of standard perturba- 
tion theory to Lagrangian-based quantum field models, the (f)\ model is superrenor- 
malizable when d = 2, 3, renormalizable when d — 4, and nonrenormalizable when 
d > 5. As seen above, the <f>\ and (p^ models have been successfully constructed 
with many properties of physical relevance verified. However, the study of (f)\ has 
made clear that there are further important subtleties to be understood in the 
construction of renormalizable quantum field models. Some of these are discussed 
in this section. 

It is a widespread view that quantum field models which are not asymptotically 
free (see below), such as 0| with positive coupling constant, are not mathematically 
consistent. A close study of (p\ shows that this view must be nuanced. Another, 
closely related point is that perturbation theory is the primary tool used by field 
theorists to make predictions that may be checked against laboratory measure- 
ments. However, it is widely believed that perturbation series in QFT are diver- 
gent, at least in models involving a bosonic field (this has been proven for P(4>)2 
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models by Jaffe [117] and for the (j>\ model by de Calan and Rivasseau [40], but the 
initial suspicion that this should be true goes back to a simple heuristic argument 
of Dyson [5S] for QED). So one can hope that perturbation theory is asymptotic 
to an exact model, or even better, since the connection between the perturbation 
theory and the exact theory is then tighter as the quantities in the exact theory 
can be uniquely recovered by the Borel procedure from the corresponding pertur- 
bation series, the perturbation series is Borel summable. The Borel summability 
of the perturbation series for a number quantities of interest has been verified in 
many of the models constructed to this date. However, an examination of (j)\ is 
revealing also in this connection] 19 ! 

Many authors have proven bounds which establish the local existence of the 
Borel transform for the standard perturbation series for <f>\ {i.e. the Borel transform 
of the series exists and is analytic in a disk centered at 0), the most recent of which 
is due to Kopper |131j (an overview and references to the earlier work may be found 
in [13 1J . as well). 

Important insights have been won in QFT through the development of renor- 
malization group techniques. These have many concrete realizations, some of which 
have been established in a mathematically rigorous manner. We briefly describe 
one of the latter type in the specific case of <$>\. Let C K be a suitable ultraviolet 
cutoff free covariance such as ( 13. 5 j) or 

C K = (-A + m2)- 1 e -(- A + m i)/ K2 , 

for which Euclidean square momenta larger than k 2 are either strongly suppressed 
or totally eliminated and for which C K — > (—A + rnfy' 1 in a suitable sense as 
K — > oo. The corresponding cutoff Euclidean measure is 



where 



e -v(x MfiCK 



V{\, A,k) = J KZl : 4 K : (x) + l -Z K 5ml : # : (x) d d x . (3.6) 

This expression involves a wave function renormalization Z K , a mass renormaliza- 
tion o~m 2 K and a coupling constant renormalization 5X K = A K — A, all of which are 
^-dependent and some or all of which divergent as k — > oo, depending on the 
value of d. These quantities can be given explicitly and are motivated by standard 
perturbation theory, as above. As seen above, when d = 2, 3 Z K may be chosen to 
be 1 and the mass and coupling constant renormalizations are polynomials in A. 
When d = 4 the three renormalization counterterms are given by power series in 



19 It should be emphasized that even when the series is Borel summable, at some order of 
the expansion the difference between the exact theory (if it exists) and the perturbation series 
prediction will get larger, not smaller, as one takes higher and higher orders into account. The 
same is true of the difference between the experimentally observed result and the perturbation 
series prediction. 
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A which are likely to be divergent. They must therefore be defined implicitly, and 
for this purpose renormalization group techniques can be applied. 

The basic idea of the renormalization group applied in this context is to break 
up any integral with respect to /^a,a,« into a sequence of integrals each over the field 
restricted to a particular range of momenta. To outline one concrete realization 
of this idea, let 5C = C K — C K > with k' < K and 5<p = (f) R — (j) K i. This split 
of the field into "high" and "low" momentum parts <p K = <p K > + 5(f) leads to the 
replacement of /za,a,k with a suitable product measure /^a,a,«;' x fJ>\,A,sc- Carrying 
out the integration over the measure on the right hand results in an effective action 
K/ (</>«') determined by setting the value of the said integration equal to 

up to a normalization factor. A computation determines the lower order in \ K 
contributions to V K >, and with a suitable choice of Z K > one finds that V K > has, 
up to "small" terms, the same form as ( I3.6p j 20 l The new renormalization terms 
Z K r, Sm 2 .,, X K i can be computed as functions of Z K , Sm 2 ., X K , resulting in the renor- 
malization group "flow equations" o In order that the "small" terms actually are 
small, it is necessary to choose k' suitably close to k and X K suitably small. Thus, 
beginning with a large value of k one must proceed in many incremental steps 
down from Kn = k to a conveniently small kq at which the final integration can 
be relatively easily estimated to provide the desired bounds. But since the limit 
K — > oo must ultimately be controlled, one would also need to have A KJV — > 
as ktv = k — ?• oo, i.e. as N — > oo. If this is so, then the model is said to be 
asymptotically free (in the ultraviolet regime). 

In cf)\ the flow equation for A n = A2™ is 

A n -i ~ A n — /?2A n — /^3A n 

with /?2 > 0, so that for small A n the preceding A n _i is smaller than A n , contrary 
to the above picture. This is a signal that the \4>\ theory is not asymptotically 
free and the renormalization group techniques cannot be applied in this case. 

However, the same flow equation shows that for A < this procedure may be 
applicable. In fact, Gawedzki and Kupiainen [76] have carried out this procedure 
to provide a rigorous construction of the (hierarchical) Euclidean A04 model for 
negative coupling constant. The standard renormalized perturbation expansion is 
asymptotic for the Schwinger functions. However, since reflection positivity is most 
likely not satisfied in this model, one does not arrive ultimately at a Minkowski 
space theory. Nonetheless, for spacetime dimensions less than 4 the corresponding 



20 These "small terms" are generally ignored in heuristic QFT; however, showing that they can 
be controlled is one of the main technical problems in the rigorous use of renormalization group 
ideas. 

21 In heuristic QFT the flow equations are usually expressed in terms of differential equations, 
e.g. the Callan-Symanzik equation, which, however, have not been given a mathematically rig- 
orous, nonperturbative basis, in general. 
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flow equation for <fi\ turns out to be consistent with the above picture also for 
positive coupling constant, and this approach yields another rigorous construction 
of the A 01 model for positive A. 

Another approach to the study of <j)^ models was motivated by Symanzik's 
insight |188] that Euclidean (p^ theory can be understood as a classical gas of 
weakly self-avoiding random paths and loops. To get some idea of this, consider 
the operator —A + m^ in its approximation as a difference matrix on the lattice 
Z d , and write it as a sum of diagonal and off-diagonal terms: 

/T 1 /- J= (2d + m 2 )I-J, 

where I is the identity matrix. The entries of the matrices J, J are indexed by the 
lattice sites in Z d , and the value of J xy is 1 if xy is a lattice bond and otherwise. 
So the covariance of the corresponding Gaussian lattice measure can be written as 



-a + mi)-,l = (r 1 / - J):i = £ ^ n+1 ( jn 



xy 
n=0 



The quantity (J n ) xy can be interpreted in this picture as a sum over all possible 
walks along the lattice which go from site y to site a; inn steps. Observe that 
(3 = ^d+ml)^ 1 < 1 and that J xy is exponentially damped as the distance between 
x and y grows to infinity. If w is a path ( a "walk" ) along the lattice using nearest 
neighbor bonds in the lattice (here viewed as a (hyper) cubic lattice) and |iu| is the 
number of lattice bonds in the walk w, then the random walk representation of 
the covariance is 

("A + m^= £ /3 H+1 (J H W, 

w.y— >x 

where the sum is over all walks from y to x. Symanzik showed how one could use 
this random walk representation for the two-point Schwinger function of the free 
field to give an expression for the Schwinger functions of the interacting model (on 
the lattice, in finite volume). This results in studying correlations of random paths, 
where the weight functions on the random paths is significantly more complicated 
than the straightforward exponential weight above |188j . A multitude of different 
random walk representations is to be found in the literature - see the monograph 
by Fernandez, Frohlich and Sokal \§7\ for a unified presentation of many of these. 

Typically, the random walk formalism is used in a finite volume, lattice approx- 
imation of the model to derive various kinds of correlation inequalities, which then 
generally carry over easily to the infinite volume limit. These correlation inequal- 
ities are utilized to establish bounds on the continuum limit, i.e. the exact model, 
and on the critical exponents of the model. Using these techniques, it has been 
definitively established by Aizenman [T] and Frohlich [72] that for d > 4 the exact 
limit theory is trivial, e.g. the exact model is Gaussian, manifest no interaction, 
has trivial scattering. On the other hand, in the case of d = 4 the situation is more 
subtle, as the limit may depend on the way the limits are taken and which details 
enter into the particular mode of construction — the reader is referred to (67] for 
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details. In many of these circumstances it has been proven that the exact theory 
is again trivial in four spacetime dimensions. This is striking, since the standard 
renormalized perturbation series for the model exists to all orders and is not triv- 
ial. Hence, the perturbation series is not asymptotic to what is apparently the 
exact theory. Moreover, the classical limit of that exact quantum theory could not 
coincide with the classical 04 field theory. 

This brings us to the question: which criteria do we use to decide when a given 
mathematical model is a/the M^ quantum field theory? For most of the models 
discussed to this point one has the reassurance that the standard (possibly renor- 
malized) perturbation series is asymptotic to some set of important quantities, 
e.g. the Schwinger functions, in the "exact model" (and in many cases the lat- 
ter quantities can be uniquely recovered from the Borel transform of the former 
series). For some of these theories one has even been able to prove that in the "ex- 
act model" the fields satisfy the (suitably interpreted) semiclassically motivated 
field equations. For further reassurance, some researchers would also like to know 
that (A) the classical limit of the "exact model" coincides with (B) the classical 
theory associated with the Lagrangian selected by whatever semiclassical reason- 
ing went into the choice of the model. This latter question has not yet received 
much attention from mathematical physicists, but Donald jl9] has established a 
result for P(<p)2 models with convex polynomial P which relates the classical limit 
of certain quantities in the quantum field model to corresponding quantities in 
the classical model. However, the relation between (A) and (B) is not a 1-1 cor- 
respondence, even in the well-behaved P{4>)i models. For example, Slade |175] 
has shown that for values a of the classical field for which the classical potential 
Uo(a) = P(a) + \m 2 a} does not equal its convex hull, the correspondence between 
(A) and (B) can break down. 

As already seen, none of these desiderata is satisfied by the (trivial) "exact 
model" for <\>\. These and other considerations have motivated Klauder to propose 
an alternative way to construct "exact" <\>\. The basic idea is not to perturb the 
free measure, introduce counterterms motivated by standard perturbation theory 
and then control a number of limits, as has been done in the work indicated above, 
but to perturb a "pseudofree" measure, introduce a different set of counterterms 
and then control certain limits to get the "exact model" [126] . This program has 
not been completed, so the reader is referred to |125U126| for more details. 

3.4 Yukawa^ models 

Theories with fermions are amenable to Euclidean methods, as well, if the in- 
teractions are quadratic in the fermions |151j . including therefore Yukawa- type 
interactions, i.e.. interactions of the form ipTip(j), T = 1,275 = —7071, for the scalar, 
respectively pseudoscalar, coupling . Initially, the Fermi fields were "integrated 
out", following Matthews and Salam, resulting in an effective action involving 
only the boson field called the Matthews-Salam determinant. Due to the neces- 
sity for a mass renormalization, the renormalized Matthews-Salam determinant is 
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given by 

det ren (l + XK(<p)) = det 3 (1 + \K(<f>)) exp 
where 



\3 f 

—tiK 3 (<P) + 5m 2 :(j) 2 : (x) d 2 x 



det n (l +A) = det [(1 + A)e^=i ^~ A)k ] , 

K = ST(j), S is the Euclidean fermion propagator and 5m 2 is a boson mass coun- 
terterm arising from second order perturbation theory. The ultraviolet cutoff index 
k has been suppressed. Seiler |170] showed that, after introducing a suitable finite 
volume cutoff into K and J : cf 2 : (x) d 2 x, the renormalized determinant has the 
necessary integrability properties with respect to the free boson measure \ic as 
k — > oo. Subsequently, Magnen and Seneor |139] and Cooper and Rosen [13] inde- 
pendently verified the Osterwalder-Schrader axioms for the Schwinger functions 
of the infinite volume limit by using cluster expansions. The Fermi fields in the 
Schwinger functions are also "integrated out" : the Schwinger function for n boson 
fields, m fermion and m anti-fermion fields is given by 

/n 
( H (f>(xi)) det [S'( yi , z k ; 0)] det ren (l + \K(<fi) dfictf) , (3.7) 



i=i 



where 



det ren (l + \K(0))dfic(<P) 



the determinant is applied to the matrix whose (i,j)ih entry is S ; (yi, Zk',(j>), i,j = 
1, . . . ,m, which is the two point Schwinger function for the fermions in the ex- 
ternal field and determined by (1 + XK)S' = S (all cutoffs are suppressed for 
transparency) . 

Renouard [157] proved that the perturbation expansion for the Schwinger func- 
tions of the limit theory is Borel summable. Balaban and Gawedzki [13] established 
the existence of a phase transition in the two-dimensional Euclidean pseudoscalar 
Yukawa model for large fermion mass. To do so, they adapted the mean field 
cluster expansion method developed by Glimm, Jaffe and Spencer [HI] to this 
situation. 

Lesniewski [137] utilized advances in CQFT developed for purely fermionic mod- 
els to reverse the initial procedure. He "integrated out" the bosonic field first, 
resulting in an effective fermionic action of the form 

2 / 9k(x -y) :%l)ip: (x) : i/yijj : (y) dxdy , 

where 

9k{x - y) 



(2ir) 2 J p 2 + m 2 + 8m 2 K (X) 

A is the coupling constant and the mass counterterm 5m 2 K (X) is given by second 
order perturbation theory. Using renormalization group techniques developed to 
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treat the Gross- Neveu model (see the next section), he reproved the existence and 
the Borel summability of the limit theory. 

By modifying the phase space cell expansion developed by Glimm and Jaffe to 
deal with the (/)% model, Magnen and Seneor [141] established the Osterwalder- 
Schrader axioms and the Borel summability of the standard perturbation theory 
for the Schwinger functions of the pseudoscalar Yukawas model and also indicated 
how to prove the corresponding results for the scalar Yukawa3 model. 

3.5 Gross— Neveu2 model 

The Lagrangian for the (massive) Gross-Neveu model is given by 

ip(x) (i($ + m)ip(x) + — (ijj(x)ip(x)) , 

where ip is a fermion field with iV components (colors) and m > is the fermion 
mass. In the original model proposed by Gross and Neveu in 1974 [96], the bare 
mass was 0, but the field "acquired a mass" by a complicated mechanism called 
dynamical symmetry breaking. This heuristic picture has been partially supported 
by results by Kopper, Magnen and Rivasseau for sufficiently large N [130], but we 
shall restrict the discussion to the massive case, for which definitive results have 
been attained. 

In two spacetime dimensions the model is renormalizable and asymptotically 
free when N > 1. The (somewhat simplified) spatially and ultraviolet cutoff action 
is given by 

r \ N N N 

J A a=l o=l o=l 

where 8m is the mass counterterm, 5( the wave function counterterm, and iV > 
1. The spinor indices and the ultraviolet cutoff subscripts are suppressed for 
transparency. 

Taking advantage of the fact that perturbation theory for fermions is much sim- 
pler than that for bosons, Gawedzki and Kupiainen [75], on one hand, and Feldman 
et alia [65], on the other, provided different proofs that the Schwinger functions 
converge, as the volume and ultraviolet cutoffs are removed, to Schwinger func- 
tions satisfying the Osterwalder-Schrader axioms for all sufficiently small values 
of the renormalized coupling constant. Moreover, the limit Schwinger functions 
are the Borel sum of their standard renormalized perturbation series [65J, i.e. the 
expansion in the renormalized coupling constant. By employing more recent ad- 
vances in the rigorous treatment of renormalization group methods, Disertori and 
Rivasseau [H] have found a significantly simpler proof of these results and have 
shown that the mentioned Borel summability is uniform in N. Moreover, they 
demonstrate that the renormalization group equations and the /3 function are well 
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defined in the fully interacting limit theory. Iagolnitzer and Magnen |11UJ stud- 
ied the Bethe-Salpeter kernel in this model and proved two-particle asymptotic 
completeness in this model. 

The Gross-Neveu model has also been studied in 3 spacetime dimensions, where 
it is neither renormalizable nor asymptotically free. Heuristically, however, it is 
asymptotically free in the limit N — > oo. After some changes of variables to 
represent the model for large N as a perturbation of its N — > oo limit and a rigorous 
summation of the simplest, most divergent coupling constant contributions that 
changes the nonrenormalizable model into a renormalizable one, de Calan, Faria de 
Veiga, Magnen and Seneor have shown that for sufficiently large N the Schwinger 
functions of GN 3 exist after the ultraviolet and volume cutoffs are removed, though 
most of the Osterwalder-Schrader axioms have not yet been verified (cf. [39]). It 
is striking that the exact Schwinger functions exist but the standard perturbation 
series in this model does not, since it is not renormalizable. 

3.6 Sine— Gordon2 and Thirring models 

In this section we discuss a set of models related to the sine-Gordon model, which 
is a model of a real, scalar massive or massless bosonic field <fi in two spacetime 
dimensions with (spatially cutoff) action 

Va = C, \ '■ cos(a0 + 9) : (x) d 2 x , 
Jk 

where a,( el and 9 £ [0,2tt). A number of different parametrizations of this 
action are to be found in the literature. 

The ultraviolet divergences of this model depend on the size of |a|. If a £ 
( — 2-y/7r, 2y/n), the model is superrenormalizable and already the measure 



MA 



eVc 



/ e^W duett) 



is well defined. In this range for a the model is amenable to the techniques 
previously developed to construct P{4>)2 models and models with exponential in- 
teraction in two spacetime dimensions. Frohlich and Seiler [73J proved that for 
sufficiently small (Cl/^o the standard perturbation series expansion in ( for the 
Schwinger functions actually converges in the infinite volume limit, the Schwinger 
functions satisfy the Osterwalder-Schrader axioms and the mass hyperboloid is 
isolated. The Haag-Ruelle scattering theory is applicable, and the S-matrix is 
nontrivial. Moreover, the perturbation series in ( is asymptotic to the scattering 
amplitudes. 

On the other hand, Park |156] showed that for the massless sine-Gordon model 
the infinite volume limits of the expectations (with respect to the corresponding 
finite volume measures) of products of (smeared) fields of the form : cos e(0 + 9) : 
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(x), : sin e(4> + 0) : (x) and a- V0(x) (a a constant vector) satisfy the Osterwalder- 
Schrader axioms (with the possible exception of clustering) for all parameter values 
in the indicated ranges by using suitable correlation inequalities. 

When a 2 £ [4ir, 8ir) the model is superrenormalizable, but the number of renor- 
malization counterterms increases to infinity as a increases to 8ir 2 . Nicolo, Renn 
and Steinmann [148 J proved that in a finite volume the massive sine-Gordon model 
is ultraviolet stable for a in this range. In the same range Dimock and Hurd use 
renormalization group methods to prove bounds on the finite volume Schwinger 
functions for both the massless and massive cases which are uniform in A in the 
massive model. Moreover, they show that in a finite volume the Schwinger func- 
tions are analytic in £ at 0. However, the full control of the infinite volume limit 
has not yet been attained. 

When a 2 = 8ir the model is heuristically renormalizable but not superrenor- 
malizable, and when a 2 > 8tt it is nonrenormalizable. Nicolo and Perfetti |149] 
have proven that for a 2 = 8tt the model is indeed perturbatively renormalizable. 
And for a 2 > 8ir Dimock and Hurd [47] have shown that the model is asymptoti- 
cally free in the infrared. However other rigorous results in these cases are not yet 
available. 

One of the first manifestations of the equivalence between apparently distinct 
quantum field models is the phenomenon of bosonization in two spacetime dimen- 
sions. A simple example is the equivalence between free massless fermionic fields 
and free massless bosonic fields given by the identifications 

^(1 + a 75 )V ++ c : e iaV ^ <t> : , ^if) ++ W^fl^, 



where a = ±1, e is the standard totally antisymmetric matrix with entries ±1, 0, 
and c is a suitable constant depending on the choice of Wick product used. Cole- 
man [42] gave a heuristic argument suggesting a similar equivalence between the 
massive Thirring model with Lagrangian 

where Z, Z\ are renormalization constants and j^ = ipj^ip, and the massless Sine- 
Gordon model with Lagrangian 



d^d^cf) + ( : cos(q;<; 



under the identifications 

Z^){\ + a 75 )V *+ c : e iaa * : , Z^^ +* -ci^d v <j> , 

where c, c\ are suitable constants depending on A and the ultraviolet cutoff used. 
In order for this equivalence to be valid, certain relations among the Thirring 
parameters A, r and the Sine-Gordon parameters (, a must obtain. The case 

25 



a 2 = 47r corresponds to free fermions (A = 0), and the case ( = (free bosons) 
corresponds to massless fermions (r = 0). 

Closely associated with these models is the Thirring-Schwinger model, often 
referred to as QED2 since the only manifestation of electromagnetism in one spatial 
dimension is the Coulomb force between charges. The interaction Lagrangian for 
this model is 

where ip is a massive two-component Fermi field, j^ is the conserved current defined 
above, and j° = j° + j®, with j® is a formal c-number current specifying charges 
at infinity, and V — ||x| is the one-dimensional Coulomb potential. Frohlich and 
Seiler [73] rigorously proved the equivalence between this model and the massive 
sine-Gordon model for sufficiently large mass and a 2 < 4tt by showing, along the 
lines of Coleman's original argument, that their perturbation series coincide term 
by term and actually converge. 

From the point of view of heuristic Lagrangian considerations, the Thirring 
model is renormalizable but not superrenormalizable. In point of fact, different 
versions of the (massless) "Thirring model" are extant. There is the original model 
introduced by Thirring |191j . for which Glaser [80] found an explicit "solution" for 
the fields. However, Ruijsenaars |161| showed that the n-point functions of the 
fields of this "solution" do not exist. Another concrete realization was proposed 
by Johnson [118] , which, although it is mathematically unsatisfactory, led to a 
rigorous formula by Klaiber |124] for Wightman functions which coincide with the 
n-point functions of Johnson's realization for n = 2,4. Carey, Ruijsenaars and 
Wright [H] proved that Klaiber's functions satisfy the Wightman axioms. In fact, 
they rigorously constructed the fields corresponding to these n-point functions as 
strong limits of certain natural approximating fields. 

A rigorous construction of the massive Thirring model was finally achieved 
by Benfatto, Falco and Mastropietro [19]. Writing the ultraviolet cutoff generat- 
ing functional for the Euclidean model as a Grassmanian integral, they showed 
that after proper choice of the wave function renormalization and bare mass the 
Schwinger functions (at noncoinciding points) converge as the cutoff is removed 
to Schwinger functions satisfying the Osterwalder-Schrader axioms. A multiscale 
renormalization group approach related to the one discussed in Section 13.31 is used 
in the proof, and so the results are valid for sufficiently small values of the coupling 
constant (there is no restriction on the mass). Curiously, they also find that in the 
massless case the resultant two-point function differs from that found by Johnson. 
Benfatto, Falco and Mastropietro [20] then went on to prove Coleman's equiva- 
lence between the massless sine-Gordon model and the massive Thirring model in 
any finite volume. 
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3.7 Local gauge quantum field theories 



Local gauge quantum field theories are conceptually and mathematically quite 
challenging, and we discuss briefly some of these challenges in this section. Al- 
though we are concerned here with CQFT, it is necessary to treat some of the 
insights gained into the nature of gauge theories by other mathematically rigorous 
means, since they shed light on the nature of these challenges and on the conse- 
quences for the mathematical framework of the models. But even the mathemati- 
cally rigorous literature on aspects of local gauge theory has become enormous, so 
only a few highlights can be discussed here. 

One of the challenges facing the rigorous construction of such models is to 
decide what constitutes the proper framework for the end result. It has become 
clear that, though the Wightman axioms are still suitable for the gauge invariant 
local fields in gauge theories, such as the electromagnetic fields F^ v , they must be 
supplemented to include extended field objects, since it is also useful to consider 
gauge invariant objects such as 

^(x)e iJc *y Mz)dz " iP(y), (3.8) 

here C xy is a suitable curve connecting the points x and y, A is a gauge potential, 
and ip a fermionic field. Various incomplete proposals have been made in this 
regard, but the most developed of these appears to be that of Frohlich, Osterwalder 
and Seiler (cf. the last chapter of |172j ). which provides conditions on expectations 
of products of Euclidean Wilson loops e^c A v( x ) dxl1 f or piecewise smooth loops 
C in spacelike hyperplanes in Minkowski space (viewed, however, in Euclidean 
space) as well as a procedure to construct the corresponding Minkowski space 
theory. Another proposal by Ashtekar, Thiemann and co-workers [TT1I190] places 
its conditions instead on the Euclidean measure of the model, instead of on specific 
classes of expectations. These conditions also afford the possibility to reconstruct 
important aspects of the real time theory from the Euclidean data. They have 
shown that a large class of Yang-Mills models in two spacetime dimensions verify 
these conditions — see Section 13.7.31 

In addition, it is convenient for various purposes to include the unobservable 
gauge potentials A^ in heuristic gauge QFT. But if they are allowed to enter into 
the theory as a quantum field in their own right, the price is high. It was recognized 
by many physicists that the introduction of gauge potentials as quantum fields 
acting on the state space is incompatible with both Lorentz covariance and Einstein 
causality. We illustrate this point with the electromagnetic field in Section 13.7.21 
The upshot is that either one must choose a gauge in which the gauge potential 
acts as an operator in a Hilbert space but is not Lorentz covariant and violates 



22 The adjective "local" modifies "gauge", not "quantum field theory". The internal symmetry 
groups in these theories are infinite dimensional, as the "gauge transformation" can depend upon 
the spacetime point at which it is being implemented. This is to be contrasted with "global gauge 
groups" , such as the global U(l) symmetry in the Yukawa models or the global Z2 symmetry in 
P(4>)2 models with even polynomial P. 
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Einstein causality, or one must choose a gauge in which the gauge potential is 
Lorentz covariant and satisfies Einstein causality but acts as an operator in a 
vector space with an indefinite inner product, i.e. there exist vectors \l/ in the state 
space such that (\1/, \1/) < 0. In such a circumstance the standard relation between 
quantum expectations and probabilities fails. One should note that perturbation 
theory is generally performed in gauges of the latter type. Additional complications 
are introduced by the (again unobservable) charge carrying Fermi fields in such 
quantum field theories (again illustrated in Section I3.7.2J) . 

A modified version of the Wightman axioms which takes into account the ne- 
cessity of employing a reasonably behaved indefinite inner product space called 
a Krein space may be found in the monograph of Strocchi [TSOj J 23 l The essen- 
tial modification is that the positivity condition on the Wightman functions (or 
the reflection positivity condition on the Schwinger functions) is replaced by a 
"HUbert space structure" condition which permits the construction of subspaces 
of the Krein space suitably associated with the Wightman functions on which the 
inner product is positive definite. 

Taking the standpoint that the Euclidean gauge potentials are to be taken 
explicitly into account as dynamical variables brings up the question of which 
measure to adopt on the space of potentials. In the geometric formulation of 
classical electrodynamics, the gauge potentials are connections on a certain fiber 
bundle with base space M. d (and the Wilson loops are (traces of) holonomies of the 
connection around closed loops). Let srf denote the set of such connections, and let 
£f denote the local gauge group acting on srf . In gauge theories it is natural to view 
each gauge equivalence class of connections as a distinct physical "path" , i. e. the 
orbit of any single potential Ae^ under the action of ^ is an equivalence class. 
The elements of each orbit are thus viewed as physically equivalent. However, 
the resultant path space — the space of all such orbits — is a nonlinear quotient 
space stf /£f ; this introduces technical difficulties. By using a gauge fixing one 
can impose a linear structure on srf ' j^ for d = 2, but in higher dimensions, the 
resultant Gribov ambiguities limit the usefulness of such gauge fixings. 

A further complication arising when one takes the gauge potentials as dynam- 
ical variables is that the measures of physical relevance in QFT are typically not 
supported on functions having nice properties from the point of view of analysis (in 
particular, the measures are not supported on £//&), so that one must introduce 
a suitable closure of &/ /& in which one can find the relevant generalized connec- 
tions on which these measures do have support. The choice of this closure is not 
unique, and many can be found in the literature. The measures are then defined 
on this closure. In the scalar and fermionic Euclidean models described above, 
the (cutoff) physical measure is defined as a suitable perturbation of a Gaussian 
measure corresponding to a free field. In CQFT this Gaussian measure replaces 



23 A modified version of the Osterwalder-Schrader axioms for the same purpose has been pro- 
posed by Jakobczyk and Strocchi J114J . They show that if the Schwinger functions satisfy these 
modified conditions, then the reconstructed Wightman functions satisfy the modified Wightman 
axioms. 
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the "measure" exp{— ml J 2 (x) dx d } du(4>)/Z to be found in heuristic treatments 
of functional integration in QFT, where v is the nonexistent Lebesgue measure 
on iS'(]R d ). For gauge models Ashtekar and Lewandowski [TD] have constructed a 
uniform Borel measure \Lq on a natural closure of srf l*S which is gauge invariant 
and strictly positive on continuous cylindrical functions. This measure replaces 
v in many rigorous treatments of functional integration for gauge theories. The 
physical measures are then obtained by suitably perturbing /z . This is illustrated 
in Section 13.7.31 The advantage of this approach is that the gauge invariance is 
explicit at every step. Any Fermi fields coupled to the gauge potentials would 
be "integrated out" along the lines of the previous work on the Yukawa models, 
resulting in a new effective action perturbing the pure gauge measure. 

The complications in defining measures on such path spaces are great, lead- 
ing many researchers to concentrate on the topological and geometric aspects of 
the problems (in so-called topological quantum field theory), largely bypassing (or 
ignoring) the analytical aspects. These interesting developments are not treated 
here. And, although also classical gravitation can be understood as a local gauge 
theory, we do not address any aspect of quantum gravitation here. However, we 
give a brief accounting of the constructive results concerning the physically most 
important quantum gauge models, even though definitive results have been at- 
tained only in spacetime dimensions less than 4 and, as indicated above, aspects 
of these models do not conform to the Wightman setting. It should be emphasized 
that it is not at all clear that gauge theories such as QED and QCD cannot be in- 
corporated into the HAK setting, since there the observable quantities are primary 
and many of the problems we have discussed above arise only after introduction 
of unobservable fields. One of the many successes of AQFT is that Doplicher, 
Haag and Roberts have shown how, starting with the (necessarily gauge invari- 
ant) observables in a global gauge theory, one can uniquely derive both the gauge 
group and the associated charge carrying fields - - cf. e.g. [9| I104] . Although a 
corresponding breakthrough for local gauge theories has not yet been achieved, 
the possibility remains open. 

3.7.1 Abelian Higgs d model 

The abelian (or U(l)) Higgs model is an interacting theory of a vector field A u (x) 
coupled in a gauge covariant manner to an A-component scalar field 4>(x) and is 
sometimes referred to as scalar electrodynamics (when N = 1). The Euclidean 
action of the model in d spacetime dimensions is 

where the field strength tensor is F^ix) = d^A v — d^A^, and the covariant deriva- 
tive of the scalar field is 

(■D M 0)i(ar) = dpfaix) - eAn(x)(Q(f))i(x) . 
29 



Q is an antisymmetric N x N matrix, and e and A are coupling constants. For 
d = 2, 3 the model is superrenormalizable. 

In a series of papers (see [31] for references), Brydges, Frohlich and Seiler studied 
this model in two spacetime dimensions by commencing with the model on a 
Euclidean lattice. For the convenience of avoiding spurious infrared divergences, 
the bare mass of the gauge field A^ is initially taken to be strictly positive, but 
ultimately the continuum limit, the ultraviolet limit and the limit in which the bare 
mass of the gauge field tends to zero are all controlled, and all the Osterwalder- 
Schrader axioms except clustering are verified for Schwinger functions involving 
gauge invariant local fields such as : \<f)\ 2 and F^ u , as well as string and loop 
observables such as 

:$(x)e% A » dx '» (j)(y) : 

and : e^ A ^ dx ^ \. Hence, there exists a corresponding Wightman theory for these 
fields, though the vacuum may not be unique. 

King [122(1123] studied the abelian Higgs model in two and three spacetime di- 
mensions. He also began with the model on a finite volume lattice and controlled 
the continuum and infinite volume limits using a combination of renormaliza- 
tion group techniques and correlation inequalities. He verified all of Osterwalder- 
Schrader axioms except ergodicity and the regularity properties which assure that 
the Wightman functions are tempered distributions. 

Balaban, Imbrie and Jaffe [15] have examined the abelian Higgs model for 
d = 2, 3 for rigorous evidence supporting the heuristic notion of the Higgs mech- 
anism, i.e. the notion that, under certain circumstances usually associated with 
spontaneous breaking of the gauge symmetry, a massless particle can acquire mass. 
A proof of the existence of a mass gap on the Hilbert space generated from the 
vacuum by application of all gauge invariant observables could be taken as such 
supporting evidence. Together with Brydges they supplied such a proof for the 
model on a lattice. Although they made much progress on the proof of the mass 
gap in the continuum limit [15], the final step in that proof never appeared in 
print. 

3.7.2 Quantum electrodynamics^ 

Despite the success of the experimental predictions made by the perturbation 
theory computations associated with QED, it is widely believed that since QED 
is not asymptotically free in four spacetime dimensions it cannot be defined as a 
mathematically rigorous theory; but this matter has not been settled. Moreover, 
beyond the ultraviolet problems inherent in the model some serious conceptual 
difficulties remain to be resolved. As a first example, though the notion of a 
gauge transformation is straightforward in classical electrodynamics, it is not at 
all clear what gauge transformations are in QEDr 4 ! and what their relation may 



24 as evidenced, for instance, by an investigation by Strocchi and Wightman |178) 
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be to the classical gauge transformations which play such an important role in the 
semi-classical reasoning commonly found in heuristic QFT. 

In addition, as shown by a number of mathematical physicists (see Strocchi's 
monograph |180j for details and references), under various sets of reasonable as- 
sumptions, the standard picture in quantum theory of a Hilbert space % serving 
as the state space for the model is inconsistent with each of the following: the 
potential field A^ is covariant, A^ satisfies Einstein causality, and Maxwell's equa- 
tions are satisfied on %. Wightman and Garding |196j have shown that for the 
free electromagnetic field a formalism due to Gupta and Bleuler is mathematically 
consistent. In particular, there exists a vector space 7i with a sesquilinear Her- 
mitian form (-,-) on which acts a unitary representation U(V + ) of the Poincare 
group satisfying the spectrum condition and operator valued tempered distribu- 
tions A^{x) and F^ u (x) = d^A u (x) — d u A fl (x), which are covariant under the action 
of U(V\) and satisfy Einstein causality. In V, is a distinguished subspace V.' such 
that (\l/, \l/) > for all $6^' and in which Maxwell's equations hold: 

($,&>F lil/ (x)V) = Q 

for all $, ^ e H'. The (pure) physical states are described by unit vectors in 
the Hilbert space given by the quotient H'/H", where %" is the subspace of %' 
consisting of vectors \l/ such that (\1/, \1/) = 0. For these the standard quantum 
probability interpretation is valid. 

Though there are a number of realizations of the Euclidean free electromag- 
netic field in spacetime dimension d, the most concise is to specify its generating 
function: 

/^expf-^a-A,)- 1 ^,^)], (3.9) 

where A d is the Laplacian on M d , (5f) u (x) = J^ df fMV (x)/dx fJl and / is a tempered 
function valued 2-form. This satisfies the hypotheses of Minlos' theorem, so there 
exists a probability measure fi on the dual space of the space of test 2-forms 
such that J exp[iF(f)] d\x equals (13. 9p . The Euclidean free electromagnetic field 
can be defined as the generalized stochastic process / h-> F(f) satisfying that 
equation. The corresponding Schwinger functions satisfy all of the Osterwalder- 
Schrader axioms (note that the gauge potentials do not appear at all). The exterior 
derivative of F is zero. Realizations of the Euclidean free electromagnetic field 
involving the gauge potentials run into the problems indicated above. 

Additional complications are introduced when one considers the electron field 
ip(x). Because the photon is massless and Gauss' Law holds, the physical electron 
field carrying the electric charge cannot be Lorentz covariant or satisfy Einstein 
causality. Moreover, in a Hilbert space Ti the electron field ip(x) cannot satisfy 
Einstein causality if Maxwell's equations with a source hold in TL. Furthermore, 
in fully interacting QED the physical state space will be complicated not only by 
the considerations mentioned above and by the superselection rule associated with 
the global electric charge operator, resulting in countably infinitely many distinct 
charge sectors, but also by the uncountably infinitely many sectors associated with 
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distinct asymptotic configurations of soft photons. (For a profound investigation 
of these matters, the reader is referred to a paper of Buchholz [23] •) In heuristic 
QED this abundance of superselection sectors is ignored, and one selects more or 
less arbitrarily a subset of physical states to work with, typically discarding in this 
manner the charged states with the best possible localization properties. 

The perturbation theory for QED is on solid mathematical footing. QED is 
renormalizable in four spacetime dimensions (see, in particular, the proof in [66J), 
so the perturbation series is well defined to all orders. Indeed, it has been proven 
by Feldman and co-authors [66] that the standard perturbation series for the 
Schwinger functions in QED4 is locally Borel summable. It is not known if there 
exists an exact model to whose Schwinger functions this series is asymptotic. So 
we turn now to lower dimensions. 

Weingarten and Challifour |193] formulated QED2 on a finite lattice and showed 
that a Salam-Matthews formula for Schwinger functions analogous to (13. 7\\ holds, 
where the reference measure is that of the lattice gauge theory discussed in Sec- 
tion [3TTT3] with gauge group G = U(l) (which gives a lattice approximation to the 
action of free Euclidean electromagnetic fields expressed in terms of the electro- 
magnetic potentials), the product of scalar fields is replaced by a polynomial of 
electromagnetic potentials and the renormalized determinant is replaced by a sim- 
ilar expression. They showed |193[ll94j that in the limit as the lattice spacing goes 
to and the volume cutoff is removed, their expressions for the Schwinger func- 
tions involving arbitrary products of the gauge potentials and products of pairs of 
fermionic and antifermonic fields have well defined limits. However, none of the 
Osterwalder-Schrader axioms were addressed. Seiler |172] has given a sketch of 
a cluster expansion whose convergence would verify all the Osterwalder-Schrader 
axioms for the Schwinger functions involving gauge invariant local fields (such as 
the electromagnetic field and the fermion current) plus a mass gap for the fermion 
(electron) for sufficiently small ratio of the electric charge to the fermion mass, but 
the details have not appeared in print. 

3.7.3 Yang-Mills d models 

There are many more or less equivalent formulations of the classical Yang-Mills 
theory. The most concise, which requires familiarity with some basic concepts 
of differential geometry, is the following: Let G be a compact group, A be a 
connection on a G-bundle over M. d (A is a 1-form taking values in the Lie algebra 
(25 of G), and F = dA + A A A be the corresponding curvature. The (classical) 
action for the pure Yang-Mills^ model is 

1 

4? 

where Tr denotes an invariant quadratic form on (5 and g plays the role of a 
coupling constant. Introducing coordinates on M. d and a basis in (25, one has 

F - ft A — ft A — r bc A , A 
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where c„ are the structure constants of (5. In these terms the action is 



W Jri *-? 



A 9 

If G = U(l), then the structure constants are and one recovers the action of 
pure electromagnetism. In Yang-Mills theories G is understood to be nonabelian. 
For example, in QCD one takes G = SU(3) and in the electroweak theory of 
Glashow, Salam and Weinberg one takes G = SU(2) x U(l) (both of these are 
components of the SM). In QCD and the SM there are, in addition, matter fields 
which are coupled to the gauge potentials A and provide a further contribution to 
the action. These are not discussed here, with one exception, since the rigorous 
results in these cases are minimal at this point in time. However, there are results 
of interest concerning pure Yang-Mills models, and we turn to those next. 

For spacetime dimension d > 2 one is faced again with ultraviolet problems 
that go beyond Wick ordering. And since almost all means of regularization (i.e. 
suppressing high enery-momentum values) destroy gauge invariance, most rigorous 
studies of (Euclidean) Yang-Mills in higher dimensions commence with the theory 
on a lattice. There are various versions of this, but we will consider only one, due 
essentially to Wilson |197j . 

Let A be a hypercubical lattice in M. d with lattice spacing a. Ordered pairs 
xy of nearest-neighbor lattice points x, y G A are called bonds, closed loops uvxy 
consisting of the obvious four bonds are called plaquettes. Let G be a compact Lie 
group, let x be a character on G, and let g. be a map from the bonds in A into G 
such that g xy = g~£. The collection of such maps is called the field configuration 
space. To each plaquette P = uvxy in A corresponds the conjugacy class gp of 
the element g U v9vxgxyg y u and the quantity Ap = \{x(.9p) + x(9p))- The Wilson 
action 25 ! on the lattice A is defined in terms of these quantities: 



A A 



-iE* 



w _ 
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where the sum runs over all plaquettes P in A. A probability measure on the field 
configuration space is given by 



-K 



_ 11 d 9xy 

xy 



Je- A A ]\dg. 



where the product is taken over all bonds xy in the lattice and dg xy is the Haar 
measure on G for each such bond. Gauge transformations are given by g xy (->■ 
'Jx9xy1y 1 , where 7 : A — > G. Both the action and the measure are gauge invariant. 
As shown by Osterwalder and Seiler |154j . the expectations with respect to /xa of 
gauge invariant functions on the field configuration space satisfy a natural analogue 



25 Another, less commonly used action for lattice gauge theories is the Villain action, which is 
an approximation to the Wilson action. Its description is lengthier and is not presented here. 
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of the reflection positivity condition in the Osterwalder-Schrader axioms, so that 
one can construct a corresponding Hilbert space and a positive self-adjoint transfer 
matrix (also bosons and fermions can be naturally incorporated into this setting 
|154j ). Much rigorous work has treated the infinite volume limit of such theories, 
but for the purposes of quantum field theory one also must control the limit as 
a — > 0. For this, there are significantly fewer results. 

Gross [97] showed that as the lattice spacing converges to the U(l) 3 lattice 
gauge model with the Villain action converges to the free Euclidean electromag- 
netic field in the sense of convergence of the characteristic functions of the field 
variables F^ u . He also showed that if the Wilson action is adopted, then the char- 
acteristic function of the "lattice current" converges to the characteristic function 
of the current of the free Euclidean electromagnetic field as the lattice spacing 
goes to 0. Driver [53] proved that similar results hold for the U(l) 4 lattice gauge 
model. 

In two spacetime dimensions, when the complete axial gauge is taken the Eu- 
clidean Yang-Mills action simplifies sufficiently to define a Gaussian measure for 
the gauge field. This greatly reduces the difficulties in the task of construction 
and affords the possibility of avoiding lattice approximations (though at the cost 
of giving up explicit gauge invariance) , and a number of different approaches have 
been developed to that end. Only a few highlights can be mentioned here. Gross, 
King and Sengupta [98] constructed the Euclidean U(iV) pure Yang-Mills model 
in the axial gauge in two spacetime dimensions using stochastic differential equa- 
tions (c/. Section |5]). They show that the differential equation controlling parallel 
transport along a smooth curve in M? with respect to a typical gauge potential 
can be interpreted as a stochastic differential equation. Without cutoffs they 
construct the Schwinger functions for the (nonoverlapping) Wilson loops, provide 
closed expressions for these Schwinger functions and verify their Euclidean invari- 
ance. Driver [53] extended their results to include expectations of more general 
functions of parallel transport along a finite "admissible collection" of curves and 
further proved that the continuum model is the limit of the (gauge fixed) lattice 
approximations (for both the Villain and the Wilson actions) as the lattice spacing 
goes to 0. Anshelevich and Sengupta [6] have shown that the limit N — > oo (with 
N/g 2 held fixed) of these Euclidean U(N) pure YM2 models exists and is described 
by a concrete free stochastic process. 

These results have been extended in many respects by Ashtekar and co-authors 
[TT] . Studying Euclidean pure YM 2 models with G equal to SU(iV) or U(l), they 
take the reference measure /io on a suitable closure of £//& discussed above, per- 
turb it by exp{— A^} (see below), normalize to obtain a probability measure, and 
explicitly compute the expectations of products of Wilson loops. They then show 
that the resulting expressions have a well defined limit as the volume cutoff is 
removed and the lattice spacing goes to Oo They therefore have rigorously es- 
tablished a closed expression for the expectations of products of generic Wilson 



26 Fleischhack |68j has repaired some technical problems in their work, recovering their results. 
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loops for YM 2 . In addition, they explicitly construct the corresponding real time 
model and establish the equivalence between the Euclidean and real time formu- 
lations. Klimek and Kondracki |128] have constructed a measure without cutoffs 
for Euclidean YM 2 with G = SU(2) coupled to a massive fermion, but they did 
not verify any of the usual axioms. 

There has been progress in the rigorous analysis of YM 4 , primarily for G = 
SU(N), though in the few places where the details of the group make a difference 
in the proof of the estimates G = SU(2) is usually taken. In a highly technical ar- 
gument incorporating renormalization group transformations and stretching over 
several lengthy papers, Balaban (see [IB] for references) succeeded in proving the 
ultraviolet stability of the model on a lattice of arbitrary bond length. This should 
be the core of a proof that the ultraviolet limit (here in the guise of the lattice 
spacing going to 0) of the expectations of products of Wilson loops exists. How- 
ever, the series of papers came to an end before exact Schwinger functions were 
constructed and their properties were verified. Unfortunately, the same is true of 
a series of papers by Federbush on the same model (see [62J for references). 

The problem of constructing a pure Yang-Mills model in four spacetime dimen- 
sions remains open. In fact, the Clay Mathematics Institute has offered one million 
dollars to anyone who succeeds in showing that for any compact simple group G 
the corresponding pure Yang-Mills model in four spacetime dimensions exists, 
satisfies conditions at least as strong as the Wightman or Osterwalder-Schrader 
axioms, is nontrivial and manifests a mass gap. 



4 Functional Integral Constructions — Minkowski 

Real time functional integrals are technically more difficult to work with than the 
corresponding Euclidean integrals, since they are oscillatory in nature and do not 
benefit from exponential suppression of large ranges of field "values", as do the 
Euclidean integrals. Although much progress has been made in the control of the 
functional integrals arising in nonrelativistic quantum mechanics, the situation for 
the real time functional integrals that are relevant to QFT is much less developed. 
The reader is referred to the monograph of Klauder |125j for an accounting of 
both. There is not yet sufficient progress to treat interacting relativistic quantum 
fields rigorously in this framework. 



5 Constructions Using Stochastic Differential Equa- 
tions 

The strong links between Euclidean QFT and probability theory have been sug- 
gested in Section [3j but we discuss another such link in this section. If B t is a 
cylindric version of 5' (M d ~ ^-valued Brownian motion, so that for each / e iS(M d_1 ) 
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b{ = B t (f) is a version of one dimensional Brownian motion, then the linear Ito 
stochastic differential equation (SDE) 



d$ = v/-A d _! + 1 £ t ° dt + dB t (5.1) 

has a stationary Gaussian solution £° with mean and expectation which coincides 
with that of the free Euclidean scalar Bose field with mass 1: 



££o (x)£ t °(y) = / 0(0, x)0(t, y) d m 

J5'(R d ) 



These stochastic processes can therefore be identified in the indicated sense. An- 
other such relation of note is provided by considering a Gaussian white noise r\ 
on S'(M. d ), i.e. rj is a 5' (M. d )- valued random variable distributed according to a 
Gaussian measure v with generating function 



/, 



e^>di/fa)=exp(--||/||*) 



For A G (0, |], the solution of the stochastic partial (pseudo-) differential equation 
(SPDE) 

(-A + 1) A A = 7/ 

is given by the stochastic convolution integral 

0A = (-A+1)- A *7/. 

4>\ is a Gaussian stochastic process of mean and covariance 

E0 A (x)0 A ( 2 /) = (-A + l)- 2A (x-y), 

so that 0i/2 can be identified with the free Euclidean scalar Bose field with mass 
1. 

To get beyond Gaussian processes (and thus free quantum fields), essentially 
two approaches have been introduced. One is to perturb the linear drift term in 
(15. ip with a nonlinear term. In this case the typical sample paths are distribu- 
tions instead of functions. The second is to change the reference Gaussian process 
(Brownian motion or white noise) into a suitable non-Gaussian process, for exam- 
ple to perturb the white noise with Poisson noise. In both one typically finds a 
stochastic process which is a weak solution of the SDE one has set up; then the 
equilibrium measure corresponding to this process is the measure sought in the 
Euclidean QFT program. 

These methods or their variants have been employed by many authors for the 
purpose of constructing quantum field models. As an example, we summarize 
work of Jona-Lasinio and Mitter |119j . Let A be a square region in R 2 , Aa be the 
Laplace operator with Dirichlet boundary conditions on A and C = (— Aa + 1) _1 - 
For e E (0, ^) consider the nonlinear SDE 

dit = -\{C- e ^ t + \C l - £ : $ :) dt + dW t (5.2) 
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with some specified initial condition t j 27 l where W t is a Brownian motion with mean 
and covariance 

EW t (f)W s (g) = (f,C x - e g) min{£, S }. 

Jona-Lasinio and Mitter construct a Markov process t which is a weak solution 
of (15.2J) . is ergodic and mixing, and satisfies 



limE,0 i (/ 1 )--.0t(/n)= / 0(/l)..-0(/n)^ A (0) 



t— >oo 



where /ia is the measure given in ( 13. ip with P(0) = 4 and C the covariance men- 
tioned directly above. Hence, they have constructed, in this sense, the 02 model 
in a finite volume. Similar constructions have been made of the P{4>)2 models, the 
sine-Gordon model (in the range a 2 < 47r) and the exponential interaction model 
in two spacetime dimensions, all with a volume cutoff (see e.g. |198j ). 

The process <f>t constructed in this manner is A-dependent. Of particular interest 
is the fact that Borkar, Chari and Mitter [25] showed that one can take the infinite 
volume limit also on the left hand side of the above equality (the processes indexed 
by A converge in a suitable sense as A /* M 2 ). Since it is already known that the 
limit of the right hand side as A /* M?) yields the Schwinger functions of the fully 
interacting §\ model, one sees that this approach also succeeds in constructing 
the exact §\ model. The methods of Borkar, Chari and Mitter should also be 
applicable to general P(4>)2 models. 

Some effort has been expended to understand renormalization from the point 
of view of stochastic quantization. However, this has not yet resulted in a con- 
struction of a quantum field model requiring an ultraviolet renormalization going 
beyond Wick ordering. 

Albeverio and co-workers have constructed a class of Euclidean random fields 
via convolution from generalized white noise and have shown that the Schwinger 
functions of these models can be analytically continued to real time "Wightman 
functions" (c/. [I] for details and references). All of the Wightman axioms hold for 
these models except possibly the positivity conditional on the family of Wightman 
functions that allows one to construct a Hilbert space in which the field opera- 
tors act. Albeverio, Gottschalk and Wu [4] have shown that their 'Wightman 
functions" satisfy a weakened condition due to Morchio and Strocchi (cf. [180]) 
motivated by the study of local gauge theories (see Section [37f|) . This condition 
assures that one can construct a Krein space, instead of a Hilbert space, in which 
the special subset T> mentioned in Section 1 is dense. In these models one may 
choose d as one likes, so this approach results in QFT models in four spacetime 
dimensions represented in Krein spaces, instead of Hilbert spaces. The scattering 
behavior in such models has also been studied by Albeverio and Gottschalk [5], 
and some of the models manifest nontrivial scattering. However, it is not clear at 
this point whether these models include any of the models of interest in elementary 
particle physics. 



27 The solution and its laws will therefore depend on the choice of i. 

28 In some instances it has been proven that the positivity condition is, in fact, violated. 
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6 Algebraic Constructions II 

Relatively recently, new insights and tools attained in AQFT have led to an im- 
portant renaissance of algebraic real time constructions of quantum field models, 
resulting in new techniques and the construction of models which either cannot be 
constructed by other known techniques or can only be so constructed with a pro- 
hibitive amount of effort. In these constructions one is not guided by Lagrangian 
QFT but rather uses other input to arrive at the models. 

For many of these constructions a particular collection of spacetime regions is 
distinguished — the so-called wedges. After choosing an inertial frame of reference 
in Minkowski space, one defines the right wedge to be the set 
Wk = {x = (xo,xi,X2,xs) G R 4 | X\ > \x \} and the set of wedges to be W = 
{XWr I A G V + }. The set of wedges is independent of the choice of reference 
frame; only which wedge is designated the right wedge is frame-dependent. Let #w 
denote the reflection on Minkowski space about the "edge" of the wedge W. The 
set {9 W | W G W} generates the proper Poincare group V + . These regions have 
been of interest since an important insight won by Bisognano and Wichmann [21] , 
namely that for any net {A(0)}oen °f von Neumann algebras locally associated 
to a Wightman theory (0, H, U, f2) in a natural manneii 29 ! the abstract modular 
objects Aw, Jw associated with the pair (A(VV),Q) by Tomita-Takesaki theory 
[26 |I185[[189] have, in fact, physical meaning. They showed that one has 

J Wr = 6 C/(i?x(7r)) , A Wfl = U{v R {2nt)) , (6.1) 

where 9 is the PCT-operator associated to the Wightman field, vn(t), t G R, is 
the one-parameter subgroup of boosts leaving the wedge Wr invariant, and Ri(ir) 
is the rotation through the angle n about the 1-axis, which is perpendicular to 
the edge of Wr (similar relations are valid for the modular objects of any wedge). 
Hence, 

J Wr A(0) J Wr = A(9 Wr O) , A% r A(0)A^ r = A(v R (2irt)0) , 

for all O. Note for later use that A Wfl = U(v R (2irt)) = e i27TtKl entails 

A{£ = e^ . (6.2) 

Though some of the significance of their insight will become clear below, we refer 
the reader to |186] for further development and references. 

A fairly general strategy [TFJ [36] to construct a quantum field theory alge- 
braically is to construct a unitary representation U(V + ) (the representation would 
be fixed by an analysis of particle masses, types and multiplicities in scattering 



29 Roughly speaking, the bounded functions of the held operators (/>(/) smeared with test 
functions / having support contained in the spacetime region O are contained in the algebra 
A(<D). 

38 



experiments) satisfying the spectrum condition_j and for a fixed wedge Wo exhibit 
an algebra (5 which satisfies the consistency conditions: 

(a) U(X) (5 U(X)- 1 C 0, whenever AW C W for A G V\. 

(b) U(X') eUiX')- 1 C <3', whenever A'W C W for A' G V\, where a primed 
algebra denotes the commutant of the algebra and a primed wedge denotes the 
complementary wedge, i.e. the maximal wedge which is spacelike separated from 
the given wedge. 

Then setting A{W) = U(X) <5U{Xy\ where A G V\ is such that W = XW 
for any W G W, one obtains an "algebra of observables" for each wedge region. 
For any other (convex) causally complete spacetime region O one can then define 
A{0) as the intersection of all wedge algebras ^4.(>V) such that O C W. The 
resulting net {A(0)}oen satisfies the conditions of isotony, relativistic covariance 
and Einstein causality. Conversely, any asymptotically complete quantum field 
theory fixes an algebra & satisfying the consistency conditions. One may therefore 
view such pairs (U, 0) as germs of quantum field models. However, at present 
one does not have a general dynamical principle by which the algebras & may be 
selected, given the representation U. In the following, we briefly outline what may 
be seen as various ways of arriving at such pairs. 

6.1 Modular localization 

Though there have been many constructions of free quantum fields using various 
techniques, Bisognano and Wichmann's breakthrough motivated Brunetti, Guido 
and Longo [29] to find what may be viewed as an intrinsic construction which we 
discuss next. We restrict our discussion to d — 4; the case d — 3 is more complex, 
as anyons are admitted |145j . One of the manifold aspects of the subtle notion 
of "particle" in QFT arises from Wigner's famous classification of the irreducible 
representations of the (covering group of the) Poincare group |195j . in which each 
such representation is uniquely determined (up to a natural equivalence) by two 
numbers that can be interpreted as the mass and the spin of the particle (if the mass 
is strictly positive; when the mass is zero, there is the additional case of "continuous 
spin"). So presented with the task to construct a quantum field describing free 
particles of given mass and spin, one begins with the corresponding irreducible 
representation U\(V+) of the proper Poincare group on a Hilbert space Hi. From 
this data one constructs a representation U{V+) on H, the bosonic Fock space with 
one-particle space Hi- Writing the representation of the boost subgroup leaving 
Wr invariant in terms of its self-adjoint generator, Ui(vji(t)) = e' ltKl , one defines, 
in light of 06.11) and (16. 2p above, 

A]ll = e« K \ J Wr = Ui(9 Wr ) 



30 Indeed, due to the fact that these wedge algebras are generally type III factors and therefore 
unitarily equivalent, it would suffice to do so on Fock space. 
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Motivated by Tomita-Takesaki theory one then defines 

Sw R = Ui(9w R )e n * = Jw fl A Wfl 
and identifies the corresponding real subspace of invariant vectors in Tti 

K-Wr — {/ e D(Sw R ) | Sw R f = /} 
(similarly for all W G W). On the exponential (or coherent) vectors 

vn! 
which are total in "%, one defines unitary operators V(f), f G 7/i, by 

V(f)e° = e-^ 2 e^ f , f EHi, 

V(f)V(g) = e-^^V(f + g),f,geH 1 . 

Now let A(yVn) denote the von Neumann algebra generated by {V(f) \ f G Ky^ R }. 
Brunetti, Guido and Longo showed that if Ui(V+) satisfies the spectrum condition, 
then (U(V\_),A(Wr)) satisfies the consistency conditions and therefore determines 
a local, Poincare covariant net {A(0)}. In the case of finite spin, the resulting 
net coincides with the net of local algebras associated with the corresponding free 
field; in the special case of massless, "infinite (or continuous) spin" representations 
U\, this provides the first construction of a quantum field model covariant under 
an infinite spin representation. Indeed, there is no quantum field of the standard 
type associated with the net in this case, since the Fock vacuum vector Q = 
e° is not cyclic for algebras A(0) when O is bounded; instead, these nets are 
generated by string-localized fields [145J. However, for any choice of positive energy 
representation XJ\ the vector Q is cyclic for A(W), for any wedge W, and for A(C), 
for any spacelike cone C. In the case that U\(J > +) is irreducible with finite spin, Q 
is also cyclic for A(0), for any double cone O. 

6.2 Models with nontrivial factorizing S— matrices 

A classic problem of obvious physical importance is the so-called inverse scat- 
tering problem, which in the context of QFT is: given the scattering matrix S 
(determined, in principle, by the measured scattering data in an experiment), 
does there exist a quantum field model whose S-matrix is the specified operator 
SI As this is an enormously difficult problem, in order to make the problem more 
manageable workers in the field have focussed their attention on the special case 
of factorizing S-matrices in two spacetime dimensions; this is the relatively sim- 
ple situation where all scattering processes reduce to (suitable combinations of) 
two-body scattering, so that specification of the two-body scattering amplitude 
completely determines the S-matrix. There is, therefore, no particle production 
in such models. The primary effort in this direction has been made in the context 
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of the form factor program {121(117611199] , in which local quantum field operators 
associated with the quantum field purportedly having the prescribed scattering 
behavior are expressed in terms of a certain algebra, the Zamalodchikov algebra. 
Rigorous formulas for matrix elements of local operators between scattering states 
have been obtained, but the computation of products of local operators at dif- 
ferent spacetime points is not under mathematical control, because infinite sums 
over intermediate states are involved. Hence, the Wightman axioms have not been 
verified in such models, apart from simple cases. 



However, Schroer [167(1168] realized that certain field operators in the Zamalod- 
chikov algebra can be interpreted as being localized in wedges. These wedge- 
localized but nonlocal quantum field operators that create covariant one-particle 
states out of the vacuum are now called polarization-free generators. It was subse- 
quently shown by Borchers, Buchholz and Schroer [21] that in more than two space- 
time dimensions the existence of (tempered) polarization-free generators defined 
on a translation- invariant, common dense domain (tacitly assumed in [167(1168] ) 
entails the triviality of the associated S-matrix. However, for constructing quan- 
tum field models with nontrivial scattering in two dimensional Minkowski space 
the idea is quite fruitful. 

In the special case of these factorizing S-matrix models the S-matrix is deter- 
mined by a single function £2 through the relation 

(s*) n (e 1 ,...,e n ) = [ J] s 2 (\e l -e k \)]* n (e 1 ,...,e n ), 

l<l<k<n 

which defines its action on a general n-particle wave function in the Hilbert space 
described below, where 9{, i — 1, . . . , n, are the rapidities of the scattered particles. 
Lechner [133H135J considers a large class of two body scattering functions S 2 (9) 
satisfying conditions arising from the requirements of unitarity, crossing symmetry 
and hermitian analyticity of the corresponding S-matrix and uses them to define 
a concrete representation of the Zamalodchikov algebra. With "Hi = L 2 (M), let "H 
be the ^-symmetrized Fock space, where the n-particle wave functions satisfy 

\&n(#i, . . . , 9i + i, 9 i: . . . , 9 n ) = S 2 (9i — 9 i+ i)^ n (8i, ...,#», i+1 , ...,9 n ). 

For the special, and admissible, scattering functions S 2 = 1, —1, "H is the standard 
boson, respectively fermion, Fock space. % admits a canonical, strongly continu- 
ous unitary representation U(V + ) of the Poincare group satisfying the spectrum 
condition. On H. also act creation, resp. annihilation, operators Z\ Z, satisfying 
the Fadeev-Zamolodchikov relations: 

Z\9)Z\9') = S 2 (9 - 9')Z\9')Z\9) 

(similarly for Z) and 

Z{9)Z\9') = S 2 {9' - 9)Z\9')Z{9) + 6(9 - &) ■ 1 
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In terms of these one defines a quantum field operator 

0(/) = Zt(/ + ) + Z(/_),/e<S(R 2 ), 

where 

f±(0) = ±-Jf(x)e ± ^ e >dx 

and p(6) = m(cosh8, sinh#), where m > is the mass. 

Though the field (j>(x) is covariant under the action of U(V]_), is densely defined, 
and is a distributional solution of the Klein-Gordon equation of mass m, it does not 
satisfy Einstein causality unless S 2 = 1, in which case 4>(x) is the usual free scalar 
Bose field. However, as observed by Schroer, when smeared with test functions 
having support in a wedge, they are polarization-free generators. 

In two dimensional Minkowski space the set of wedges decomposes into two com- 
ponents, one is the set of all translates of W R and the other is the set of all trans- 
lates of Wl = W^j. Defining A(W R ) to be the von Neumann algebra generated by 
{eW> I supp(/) C W R }, A(W L ) = A(W R )', and A(W) = U(x)A(W R )U(x)- 1 if 
W = W R + x (similarly for W = W L + x), then {A(W)} We yy satisfies the HAK 
axioms for the restricted collection TZ = W. The Fock vacuum vector Q is cyclic 
for the wedge algebras and the corresponding modular objects coincide with those 
found in the Bisognano-Wichmann setting. 

In two spacetime dimensions double cones can be defined simply as the inter- 
section of two suitable wedges: O = Wi fl W 2 with W2 C Wi. Lechner showed 
that if one defines the corresponding observable algebra to be 

A(0) = A(Wi) n A(m)' (= A(Wi) n A(W 2 )) , 

one finds that for a large class of S2 (i. e. for a large class of models of the type we are 
discussing), fi is cyclic and separating for all double cone algebras and the double 
cone algebras satisfy Einstein causality. Moreover, the Haag-Ruelle scattering 
theory can then be applied to yield a scattering theory which is asymptotically 
complete and whose S-matrix coincides with the initially prescribed S-matrix 
[135j . These models therefore constitute a complete and satisfying solution to the 
inverse scattering problem for the stated class of S-matrices. 

Note that the construction is implemented using easily constructed, but nonlo- 
cal fields to obtain local wedge algebras, of which suitable relative commutants pro- 
vide algebras of observables localized in bounded spacetime regions. This sidesteps 
the usual process of first constructing local fields and then using them to obtain 
the local observable algebras. This is a significant simplification for the models in 
question, since arguments by Smirnov and Schroer |168[ll76j and examples labori- 
ously computed by McCoy et alia [144) indicate that the local fields in these models 
must be infinite power series in the (simple) nonlocal fields. By constructing the 
local algebras in the indicated manner, one is able to avoid controlling the infinite 
expansion that yields the local fields and yet still arrive at the desired quantities 
of physical relevance. 
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The special case of S 2 = — 1 was studied by Buchholz and Summers [35] for 
d > 2, and it was shown that the model is maximally nonlocal in a certain specific 
quantitative sense. For d = 2 it was shown that there are two associated local nets 
admitting an asymptotically complete scattering theory, one describing a fermion 
with trivial scattering and another describing a boson with S = (— l)^^ -1 )/ 2 , 
where N is the number operator. In higher dimensions there exist string-localized, 
respectively brane-localized, operators which mutually commute at spacelike sep- 
aration. 



6.3 Deformations 

Deformation techniques have long been used to provide quantized versions of clas- 
sical theories (see e.g. P-.92J), but attempts to deform quantum field models have 
appeared fairly recently. Their initial impetus was provided by the desire to deform 
quantum fields on Minkowski space to arrive at quantum fields on Moyal space, or 
noncommutative Minkowski space. Though most of that work has not been math- 
ematically rigorous, Grosse and Lechner [99] rigorously defined a deformation of 
the scalar massive free Bose field in Wightman's setting, which could then be in- 
terpreted either as a model on noncommutative Minkowski space or as a model 
on Minkowski space. Of particular interest is that the deformed field manifests 
nontrivial scattering (see below). Buchholz and Summers [36J then found a defor- 
mation of essentially any algebraic quantum field model which coincides with the 
deformation of Grosse and Lechner when restricted to the net associated with the 
mentioned free field. We discuss this deformation and the resulting models here. 
Though this procedure may be carried out for d > 2, we restrict our discussion to 
d = 4. 

Let ({A(0)}o£iii T~Li U, fl) satisfy the conditions stated in Section [TJ let U(x) = 
e lPx , x G R 4 , and P = J p dE(p) be the joint spectral decomposition of the 
generators of [/(M 4 ). The support of the projection- valued measure E is V + , by 
the spectrum condition. Let M be a bounded operator and, for convenience, let 
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a p (M) = e iPp Me~ iPp E B(H) ,p G 
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for fixed C, > 0, rj G R, are uniquely distinguished by the following properties [99] : 

(i) QV^cWr. 
(ii) If A = (A, x) G V\ satisfies XWr C W r , then AQA T = Q. 

(iii) If A = (A, x) G V\ satisfies XW R dW R , then AQA T = -Q. 
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The warped convolution of M is given by the formal expressior 

M Q = J dE(p) a Qp (M) , M e B{H) . (6.3) 

A detailed explication of this unusual operator valued integral is given in [37] . 
One finds that Mq is also a bounded operator on Ti, so one may further define 
the deformed algebra Aq(Wr) corresponding to the algebra A(Wr) to be the 
von Neumann algebra generated by {Aq \ A G A(Wr)}. Then (U,Aq(Wr)) 
satisfies the consistency conditions j36][37] and therefore generates, as above, a 
net satisfying the Haag-Araki-Kastler axioms. Although the original net and the 
deformed net are not isomorphic, the modular objects corresponding to the pair 
(A(VVr),Q) coincide with those corresponding to (Aq(Wr),CI) [37]. 

There are indications that the deformed algebras corresponding to bounded 
spacetime regions O may be trivial, i.e. are multiples of the identity operator, 
so one may not be able to formulate a full scattering theory for the deformed 
net. However, two-body scattering for the deformed model is well defined |36j . 
The relations between the two-body scattering states in the original and in the 
deformed theory are most transparent if one uses improper single particle states 
of sharp momentum p = (a/j> 2 + m 2 ,p), q = (\/q 2 + m 2 , q). There one has [SB] 

\ p ® Q g) in = e *W«l \p ® g) in 
\p ® Q g) out = e' l \ pQq \ \p ® q) mt . 

The scattering states in the deformed theory depend on the matrix Q through 
the choice of the wedge Wr and thus break the Lorentz symmetry in d > 2 
dimensions. This can be understood if one interprets the deformed theory as living 
on noncommutative Minkowski space, where the Lorentz symmetry is broken. 

The kernels of the elastic scattering matrices in the deformed and undeformed 
theory are related by 



out 



(p ®q q\p' Oq q') in = e ilpQqMp ' Qq ' 1 out (p ® q\p' ® qj n . 



Thus they differ, and even if the initial model has trivial scattering, the deformed 
theory does not. Hence, this deformation applied to the net of the scalar massive 
free field results in a mathematically rigorous quantum field model with nontrivial 
scattering, apparently the first such model in four spacetime dimensions. 

Dybalski and Tanimoto [57] have shown that application of the warped convo- 
lution to any chiral conformal QFT in two spacetime dimensions results in a model 
which has Lorentz invariant nontrivial scattering and is asymptotically complete, 
the first example of a massless QFT having such properties. 

These deformations were subsequently reinterpreted in the Wightman frame- 
work by Grosse and Lechner |100j in terms of a deformed product on the Borchers 



31 Strictly speaking, this is well denned for those bounded M which are smooth with respect 
to the action a x , x £ R 4 [37]. 
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algebra associated with the polynomials in the field operators mentioned in Sec- 
tion 1. And, again in that framework, Lechner [136] has found a larger class of 
deformations which includes warped deformations as a special case. In four space- 
time dimensions, the resulting deformed theories have properties similar to those 
of the warped models; however, in two dimensions he showed that the models dis- 
cussed in Section 16.21 are obtained from the free field by a deformation of this new 
type. Common to all of these deformed models so far is the absence of particle 
production. 



7 Perturbative AQFT 

As only briefly indicated in Section 1, the algebraic approach to QFT may com- 
fortably consider local algebras of observables which are not algebras of bounded 
operators acting on a Hilbert space. In this section we discuss briefly an interesting 
example of this which goes well beyond the polynomial algebras of the Wightman 
approach. This work, carried out primarily by Brunetti, Diitsch and Fredenhagen, 
explicitly incorporates perturbative QFT into the framework of AQFT. Though 
one of the primary recent concerns of these authors is to formulate perturbative 
AQFT in a manner independent of any background space-time, we shall present 
an early version of their formulation, since we are concerned here with relativistic 
QFT on Minkowski space. This will simplify the discussion somewhat. For the 
same reason, we shall not try to describe the steps taken by the authors to incor- 
porate more general interactions and shall, instead, sketch an early version of this 
approach. 

As with Borchers algebras, the connection between an abstract (complex, uni- 
tal) *-algebra A and operators acting on a Hilbert space rl (where the probabilistic 
interpretation familiar from quantum theory applies) is implemented by a repre- 
sentation, a map tc from A into the set of (densely defined) operators acting on % 
which satisfies ir(cA + B) = ck{A) + tt(B), tx(AB) = tt(A)it(B) and tt(A*) = tt(A)*, 
for all A, B G A and c G C If one finds a state u on A (a linear map u : A — > C 
such that u(A*A) > for all A G A and oj(I) = 1, where / is the unit in A), 
then such a Hilbert space "H and representation ir, called the GNS representation 
associated with u, can be canonically constructed from the data (A,u). In this 
space H, there is a unit vector Q such that u>{A) =< Q, ir(A)Q > for all A £ A. 

In perturbative AQFT (as in some forms of deformation quantization [192]) the 
basic objects are understood as formal power series. Consider the set C[[t]] of all 
complex sequences {co, ci, . . .}, where to each such sequence corresponds a formal 
power series in the formal parameter t 



oo 
C = 

n=0 



£c n f\ 



On this set addition is defined termwise and a product is defined using Cauchy's 
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formula: 

oo oo oo n 

ab=(Y,d n t n ) (J2 b n t n )=J2 (J2 a m b n-m)t n . 

n=0 n=0 n=0 m=0 

With these operations, C[[t]} is an associative, commutative ring with unit. If A 
is a complex algebra, then A[[t]] is defined as the set of all sequences with entries 
in A and, again, may be thought of as the set of all formal power series with 
coefficients in A. With addition, multiplication and scalar multiplication (with 
scalar s in C [[£]]) defined similarly, A[[t]} is an algebra over C [[£]]. 

Perturbative AQFT is based upon the causal perturbation theory developed 
primarily by Stiickelberg, Bogoliubov, Shirkov, Epstein and Glaser. This the- 
ory cannot be described here (cf. e.g. [59|), but it is a mathematically rigorous 
renormalization theory for quantum field theoretical perturbation series which pre- 
serves and utilizes in an essential manner "locality" properties, in the double sense 
commonly understood in AQFT - - localization properties coupled with Einstein 
causality. In the papers discussed below the ultraviolet problems with the formal 
series representations of the observables are handled by techniques from causal 
perturbation theory. These techniques are not described here. 

In the Fock space of the free scalar massive Bose field, consider the cutoff 
interaction Hamiltonian 



Hi{t) = - / g(t,x)A(t,x)d 3 x, 

with g an infinitely differentiable function with compact support and A a (deriva- 
tive of a) Wick polynomial. The corresponding time evolution operator from time 
— r to time r, where r > is sufficiently large that (— r, r) x M. 3 contains the 
support of g, is formally given by the Dyson series 

S(g) = 1 + 2J — / T(A(xi) ■ ■ ■ A(x n ))g(xi) ■ ■ -g{x n ) dx x --- dx n , (7.1) 



ra=l 



where the time ordered products T{A{x\) ■ ■ ■ A(x n )) are operator valued distribu- 
tions on V satisfying 

T(A( Xl ) ■ ■ ■ A{x n )) = T{A{ Xl ) ■ ■ ■ A{x k )) T{A{x k+l ) ■ ■ ■ A{x n )) 

whenever all of x^+x, ■ ■ ■ , x n do not lie in the forward lightcone of any of the points 
Xi, . . . ,Xk- Supplying the test function g with a factor (coupling constant) A, the 
local S-matrix ( 17. ip is to be understood as an element of .4[[A]], where the algebra 
A is a natural extension of the Borchers class of the free field. S(g) has an inverse 
in »4.[[A]] of the form ( 17. ip with i replaced by —i and the time ordered products 
are replaced by "antichronological" products 

T(A( Xl ) ■ ■ ■ A{x n )) = Y, ("l) |P|+n II T (Axi),i e P) , 

PeP({i,...,n}) peP 
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where V({1, . . . ,n}) is the set of all ordered partitions of {1, . . . ,n} and \P\ is 
the number of subsets in P. The "antichronological" products satisfy anticausal 
factorization. As observed by Il'in and Slavnov |112j . the local S-matrices satisfy 

S(f + 9 + h) = S(f + g)S(g)- l S(g + h) , (7.2) 

whenever the support of h is disjoint from the causal future of the support of / 
(independently of g). Further solutions of (17. 2p are obtained by introducing the 
relative S-matrices 

S^^Sigy'Sig + f), 

which also satisfy local commutation relations [S g (h), S g (f)} = whenever the 
supports of h and / are spacelike separated. 

Local algebras of observables are then introduced by defining A g (0) to be the 
*-algebra generated by the relative S-matrices S g (h) whose test functions h have 
support contained in O. If g = g' in a neighborhood of a causally closed region 
containing O, then there exists a unitary V G .4 [[A]] such that VSg^V^ 1 = S g '(h) 
for all test functions with support in O. Capitalizing upon this fact, Dutsch and 
Fredenhagen [56J construct for a given interaction Lagrangian £ (here a polynomial 
in the field or derivatives of the field) a net of such local *-algebras which satisfies 
Einstein causality and is covariant under a natural adjoint action of the usual 
unitary representation of the Poincare group defined on the initial Fock space. 

Dutsch and Fredenhagen [55] have shown that this approach can be employed 
also for QED. Because the basic quantities are local and independent of the test 
function g in the sense indicated above, no reference need be made to the adiabatic 
limit. The incompatibility in QED between gauge invariance, locality and positive 
definite inner products on the state space when the unobservable gauge potentials 
and charge carrying fields are employed is addressed by a local construction of the 
observables and of the physical Hilbert space in which the observables are faithfully 
represented (once again as formal power series of unbounded operators) . 

More recently, Brunetti, Dutsch and Fredenhagen [30] have refined this ap- 
proach to allow the treatment of low dimensional theories and non-polynomial 
interactions. In addition, they studied in this framework three of the various ap- 
proaches to renormalization group ideas that are in use among theoretical physi- 
cists and have established their mutual logical relations. They obtain an algebraic 
form of the Callan-Symanzik equation and compute the /3 function in the cj)\ and 
0g interactions in their framework, finding perfect agreement between their results 
and those found by heuristic methods. 

Although this work does not provide quantum field models satisfying the HAK 
or Wightman axioms, it falls comfortably within the framework of AQFT, since 
the primary objects are, again, nets of local *-algebras generated by observables 
which are Poincare covariant and satisfy Einstein causality and, again, the work is 
carried out with complete mathematical rigor. However, when these authors speak 
of representations in Hilbert spaces, the Hilbert spaces are vector spaces over the 
field C[[A]], not over C. So taking expectations of observables in states in this 
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approach results in a formal complex power series, not a complex number. Hence, 
in order to make the connection to experiments one must deliberately consider 
a partial sum of this series, i. e. consider the perturbation series only to a finite 
order, as is done in heuristic QFT. Since these series are not convergent, one is 
returned to the question "Is there an exact model?" 



8 Outlook 

It is evident that the efforts of the constructive quantum field theorists have been 
crowned with many successes. They have constructed superrenormalizable mod- 
els, renormalizable models and even nonrenormalizable models, as well as models 
which fall outside of that classification scheme since they apparently do not cor- 
respond to some classical Lagrangian. And they have found means to extract 
rigorously from these models physically and mathematically crucial information. 
In many of these models the HAK and the Wightman axioms have been veri- 
fied. In the models constructed to this point, the intuitions/hopes of the quantum 
field theorists have been largely confirmed. However, local gauge theories such as 
quantum electrodynamics, quantum chromodynamics and the Standard Model - 
precisely the theories whose approximations of various kinds are used in a cen- 
tral manner by elementary particle theorists and cosmologists - - remain to be 
constructed. These models present significant mathematical and conceptual chal- 
lenges to all those who are not satisfied with ad hoc and essentially instrumentalist 
computation techniques. 

Why haven't these models of greatest physical interest been constructed yet (in 
any mathematically rigorous sense which preserves the basic principles constantly 
evoked in heuristic QFT and does not satisfy itself with an approximation)? Cer- 
tainly, one can point to the practical fact that only a few dozen people have worked 
in CQFT. This should be compared with the many hundreds working in string the- 
ory and the thousands who have worked in elementary particle physics. Progress 
is necessarily slow if only a few are working on extremely difficult problemso It 
may well be that patiently proceeding along the lines indicated above and steadily 
improving the technical tools employed will ultimately yield the desired rigorous 
constructions. It may also be the case that a completely new approach is required, 
though remaining within the CQFT program as described in Section [TJ something 
whose essential novelty is analogous to the differences between the approaches in 
Section El El and EJ 

It may even be the case that, as Gurau, Magnen and Rivasseau have writ- 
ten [103J, "perhaps axiomatization of QFT might have been premature"; in other 
words, perhaps the Wightman and HAK axioms do not provide the proper mathe- 
matical framework for QED, QCD, SM, even though, as the constructive quantum 



32 Indeed, many workers in CQFT have chosen to take the methods developed for the purposes 
of CQFT and to apply them instead to problems in statistical mechanics, many body physics 
and solid state physics, where progress has been much easier. 
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field theorists have so convincingly demonstrated, that framework is quite suit- 
able for so many models of such varying types and, as the algebraic quantum field 
theorists have just as convincingly demonstrated, that framework is flexible and 
powerful when dealing with the conceptual and mathematical problems in QFT 
which go beyond mathematical existence. But it is possible that the mathemat- 
ically and conceptually essential core of a rigorous formulation of QFT that can 
include the missing models lies somewhere else. Certainly, there are presently 
many attempts to understand aspects of QFT from the perspective of mathemati- 
cal ideas which are quite unexpected when seen from the vantage point of current 
QFT and even from the vantage point of quantum theory itself, as rigorously for- 
mulated by von Neumann and many others. These speculations, as suggestive as 
some may be, are currently beyond the scope of this article. 
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